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devices.  The  generalization  is  important,  introducing  the  notion  of  incoherence  between  quantum  1/f 
fluctuations  of  the  phonon  loss  rate  in  various  regions  of  the  crystal. 

The  applications  also  include  calculation  of  quantum  1/f  noise  in  gallium  nitride.  It  is  3-10  times 
lower  in  GaN  than  in  GaAs. 
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tremendous  importance  is  the  discovery  during  this  grant  of  a  new  method  connecting  the  coherent  and 
conventional  quantum  1/f  effects.  A  mass  distribution  was  found,  which  allows  to  find  the  quantum  1/f 
noise  in  general. 

A  method  of  gate  current  suppression  in  HFET  and  a  two-dimensional  all-optical  TDM  system 
were  studied. 
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SUMMARY 


The  present  report  contains  practical  applications  of  the  Quantum  1/f  Effect,  further 
development  of  the  Quantum  1/f  Theory,  and  contributions  not  directly  based  on  quantum  1/f 
noise.  The  application  to  quartz  resonators  has  been  improved  by  the  inclusion  of  crystal  defects  in 
the  quantum  I/f  noise  calculation,  and  has  been  generalized  to  the  case  of  low-Q  and  surface 
acoustic  wave  devices,  where  Q  is  the  quality  factor.  The  generalization  is  of  great  practical  and 
theoretical  importance,  because  it  introduces  the  notion  of  incoherence  between  the  quantum  1/f 
fluctuations  of  the  phonon  loss  rate  of  various  regions  in  the  volume  of  the  crystal.  The 
applications  also  include  calculation  of  quantum  1/f  noise  in  gallium  nitride,  a  material  an  increased 
number  of  uses  and  with  a  bright  future  in  microelectronics  and  opto-electronics.  The  level  of  1/f 
noise  is  found  to  be  3-10  times  lower  in  GaN  than  in  GaAs. 

The  quantum  1/f  theory  was  reformulated  to  put  in  evidence  the  manifest  entropy 
conservation  of  the  quantum  1/f  fluctuation  process  with  the  help  of  the  Quantum  Information 
Theory  and  the  notion  of  negative  quantum  conditional  entropy.  This  reformulation  of  the  theory 
on  the  basis  of  the  new  negative  entropy  concept  explains  the  apparent  entropy  production  with  the 
help  of  a  simultaneous  production  of  negative  entropy  soft  photon  states,  A  long-standing 
conceptual  difficulty  is  eliminated  on  this  basis. 

Of  tremendous  practical  and  theoretical  importance  is  the  discovery  during  this  grant  of  a 
new  method  of  connecting  the  coherent  and  conventional  quantum  1/f  effects.  A  mass  distribution 
has  been  found,  which  allows  to  find  the  quantum  1/f  noise  in  general. 

Finally,  a  method  of  gate  current  suppression  in  HFET  with  high  dielectric  constant  gate 
insulation  and  a  two-dimensional  all-optical  time-division  multiplexing  system  were  studied. 
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I.  INTRODUCTION 


The  present  Technical  Report  contains  a  brief  presentation  of  the  main  achievements 
obtained  by  the  author  during  the  three  years  of  the  AFOSR  grant  F49620-94- 1-0342.  These 
achievements  are  both  of  practical  and  of  theoretical  nature,  the  close  interconnection  of  these  two 
aspects  being  a  main  characteristic  of  our  research. 

Our  account  will  start  with  new  [jrucilcal  results  in  the  high-tech  field,  specifically,  with  a 
necessary  generalization  of  the  author's  quantum  1/f  theory  of  frequency  fluctuations  in  the  case  of 
ultra-high-Q  quartz  resonators  to  the  important  case  of  intermediate  and  lovv-Q  resonators.  The 
first  case,  generalized  in  the  present  grant  period  to  crystals  with  lattice  defects,  is  essential  for  the 
best  resonators  at  NIST,  the  National  Institute  of  Standards  and  Technology.  They  help  define  the 
limits  of  our  notion  of  time  and  frequency,  and  we  present  this  case  in  Sec.  II.  The  second  case  is 
the  basis  of  our  industrial  applications  of  bulk  acoustic  wave  (BAW)  and  surface  acoustic  wave 
(SAW)  quartz  resonators.  It  is  given  in  Sec.  III. 

New  III-V  iTiateiials  such  as  GaN  are  of  increasing  impurlance  for  KFETs.  GaN  has  been 
used  to  produce  semiconductor  lasers  and  LEDs  over  wide  frequency  domains,  up  to  the  blue 
i*^^iOii.  ilI-\  mcitt.ricii3  ai  V  tiiisw  used  in  advanced  infrared  imaging  array’s  and  in  novel  piezoelectric 
Ill-V  semiconductor  devices.  All  these  new  materials  and  devices  are  affected  by  the  conventional 
and  often,  if  they  are  larger,  also  by  the  coherent  quantum  1/f  noise  effect.  During  the  previous 
grant  period,  the  conventional  quantum  1/f  effect  was  calculated  in  detail  as  a  function  of 
temperature  and  doping  level  for  Si  and  GaAs,  as  reponed  in  the  1994  Final  Report.  In  order  to 
help  in  the  development  of  the  new  technology,  as  shown  in  Sec.  IV.,  we  have  calculated  the 
quaiiLurti  l/f  uuise  in  GaN,  thereby  enabling  the  quantum  l/f  optimization  process  for  this  material. 

Fundamental  progress  in  understanding  the  nature  of  the  measurement  process  and  in 
eliminating  a  basic  contradiction  of  quantum  mechanics  was  recently  obtained  on  the  basis  of  the 
new  Quantum  Information-Theoretical  Approach.  This  fundamental  progress  is  based  on  the  new 
notion  of  Negative  Quantum  (von  Neumann)  Entropy  of  quantum  entangled  states.  The 
contradiction  between  the  well-known  unitary  character  of  the  evolution  of  the  state  of  a  system, 
which  precludes  any  increase  in  entropy,  and  the  actual  increase  of  the  entropy  of  the  system  as  a 
result  of  the  measurement  process  has  been  eliminated  for  the  first  time. 

The  present  author  has  applied  tiiis  fuiiuameutaily  new  concept  of  negative  entropy  for  the 
first  time  to  the  conventional  quantum  1/f  effect,  as  shown  in  Sec.  V.  As  a  result,  a  way  to 
CiiiTiinatC  tiic  long-Sttiudiiig  cwmiadiction  of  the  entropy  of  quantum  1/f  noise  appearing  from 
nowhere  was  found.  The  entropy  of  quantum  1/f  noise,  which  was  thought  to  violate  the  unitarity 
of  the  evolution  of  the  combined  system  of  electromagnetic  field  and  charged  particles  carrying  a 
current,  has  now  been  shown  by  us  to  be  fully  compensated  by  the  negative  entropy  of  the 
entangled  electromagnetic  bremsstrahlung  field.  No  other  known  phenomenon  was  so  deeply 
affectedjn  its  understanding  by  the  notion  of  Negative  Quantum  Entropy  as  the  Quantum  1/f 
niiect.  1  lie  uerivatiuii  of  the  convexiiional  quaruuiii  1/f  effect  and  of  the  Quantum  1/f  Entropy, 
without  the  use  of  any  assumptions  about  the  non-unitary  character  measurement  process,  is 
presented  in  Sec.  V. 

The  quantum  1/f  theory  provides  answers  for  any  questions  concerning  the  fundamental  1/f 
noise,  just  as  Maxwell's  equations  answer  for  electromagnetism.  Therefore,  the  quantum  1/f 
theory  must  also  provide  the  correlations  of  any  order  as  well  as  the  joint  1/f  noise  distributions 
pd  momenta  of  any  order  for  any  number  of  time  arguments.  This  ultimate,  exhaustive  and  total 
information  is  contained  in  the  characteristic  functional  of  the  process,  derived  for  subtle  case  of 
liiemiai  noise  in  Sec.  vi.  The  churacterisiic  funciiuiial  for  quaiiluiu  1/f  current  noise  was  derived 
earlier.  Finally,  we  note  that  the  fractal  dimension  of  a  1-dimensional  model  of  the  quantum  1/f 
effect  was  derived  during  the  previous  grant  period,  as  shown  in  the  1994  Report,  but  the  general 
case  has  not  been  solved  so  far  due  to  mathematical  complexity. 

A  long-standing  fundamental  problem  in  the  quantum  1/f  theory  is  the  continuous  transition 
between  the  coherent  and  conwntional  quantum  effects  for  intermediate  sizes  of  the  semiconductor 
or  metallic  sample  or  device.  The  main  problem  was  that  in  spite  of  the  obvious  relation  of  the  two 
effects  as  two  aspects  of  the  same  phenomenon,  it  was  not  clear  how  to  improve  on  the  well- 
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known  1985  interpolation  result,  or  even  how  to  start.  A  possible  practical  way  of  proceeding  was 
developed  during  this  grant  and  is  presented  in  Sec.  VII.  This  represents  a  new  way  of  bridging 
the  intermediary  region  of  device  sizes. 

Other,  less  fundamental  progress,  has  been  achieved  during  this  grant  in  the  following 
directions:  1)  Application  of  ferroelectric  and  high  e.  paramagnetic  materials  to  the  design  of  high- 
transconductance  HFETs  with  low  1/f  noise  and  low  gate  leakage  current,  as  demonstrated  in  Sec. 
VII;  2)  Discovery  of  a  two-dimensional  all-optical  parallel-to-series  TDM  and  de-multiplexing  by 
spectral-holographic  method  for  radical  noise  reduction  and  better  channel  capacity  utilization  -- 
this  one  is  also  fundamental,  but  from  a  practical  point  of  view,  as  we  see  in  Sec.  VTII; 
3)  Development  of  quantum  1/f  formulas  for  efficient  FET  and  HFET  Quantum  1/f  noise 
oplirnizalion;  4)  Development  of  quantum  l/f  formulas  for  the  three  independent  noise  sources  in 
BJTs  and  HBTs,  their  integration  in  master  formulas  and  their  use  for  efficient  device  optimization; 

5)  Development  of  methods  for  base  resistance  evaluation  from  l/f  noise  measurements; 

6)  Quantum  1/f  optimization  formulas  for  FET  and  HFET  materials;  7)  Quantum  1/f  material 
optimization  formulas  for  junction  type  devices:  junctions,  BJTs  and  HBTs;  8)  Quantum  1/f 
optimization  of  amplifiers;  9)  Quantum  1/f  optimization  of  oscillators;  and  10)  Quantum  1/f 
optimization  of  laser  clocks.  1 1)  Discovery  of  the  "thermo-electrochemical  effect".  These  points 
are  discussed  briefly  below  in  relation  to  the  corresponding  papers  mentioned  in  the  list  of 
publications  during  the  grant  in  Sec.  XI: 

1 .  The  nature  of  the  transition  between  the  coherent  and  conventional  quantum  1/f  formulas 
has  been  clarified  and  described  for  the  first  time  on  the  basis  of  an  effective  quantum  1/f  rest  mass 
distribution  of  the  electron.  This  allows  us  to  describe  in  detail  how  the  coherent  quantum  1/f 
contribution  is  effectively  reduced  to  conventional  quantum  1/f  noise  when  the  size  of  a  sample  is 
decreased.  Refer  to  papers  #7  and  3  in  the  list  of  publications  during  the  grant. 

2 .  To  provide  a  service  to  the  scientific  community  at  large,  the  nature  and  origin  of  the 
excess  heat  observed  even  in  the  "closed  cell"  variant  of  electrolytic  "cold  fusion"  experiments  has 
been  discovered.  The  effect  can  now  be  calculated  with  a  simple  formula,  and  the  "hidden 
variable"  which  caused  the  effect  to  be  absent  sometime,  was  clearly  identified.  The  quantum  1/f 
philosophy  was  applied  to  this  case.  Like  in  the  quantum  1/f  effect,  the  hidden  interaction  of  the 
system  with  the  rest  of  the  world  turned  out  to  be  the  origin  of  this  elusive  excess  heat  effect.  The 
new  effect  discovered  in  this  process  was  named  by  the  author  as  the  "thermo-electrochemical 
effect"  which  will  require  considerable  renormalization  of  many  electrochemical-thermodynamical 
data  and  tables.  It  causes  the  electrochemical  closed  system  cell  to  work  like  a  heat  pump,  even  if 
there  is  no  circulation  of  electrolyte  or  coolant  (the  system  being  closed,  self-contained  and  ideally 
insulated).  Refer  to  papers  #16  and  2  in  the  list  of  publications  during  the  grant.  The  discovery 
caused  a  partial  reconciliation  of  the  physics  and  chemistry  communities  in  the  "cold  fusion" 
debate,  a  reorientation  of  most  "cold  fusion"  efforts,  and  a  reduction  in  funding  for  those  involved 
in  this  field. 

3 .  The  quantum  1/f  effect  was  applied  to  atomic  frequency  standards,  such  as  the  Cs  and 
Rb  clocks,  yielding  a  quantum  1/f  flicker  floor  one  or  two  orders  of  magnitude  below  the  observed 
level  of  fluctuations.  This  indicates  that  the  quantum  limit  has  not  yet  been  reached  in  the  case  of 
these  clocks.  A  similar  situation  was  found  in  the  case  of  H-maser  clocks,  and  for  the  new  laser 
clocks.  All  of  these  devices  are  operating  at  fluctuation  levels  well  above  their  quantum  1/f  limit. 
Refer  to  papers  #4, 8,  and  10  in  the  list  of  publications  during  the  grant. 

4 .  A  study  of  possible  approaches  to  reducing  the  gate  current  of  HFETs  while  increasing 

their  transconductance  and  reducing  their  noise  has  shown  that  high-8  gate  insulation  provides  the 
best  solution.  This  result  was  presented  at  the  WOFE  workshop  (Tenerife,  1997),  and  has 
motivated  exploratory  research  into  the  practicality  and  optimal  realization  of  this  solution.  Refer  to 
paper  #12  in  the  list  of  publications  during  the  grant. 

5 .  A  two-dimensional  all-optical  multiplexing  and  de-multiplexing  scheme  has  been 
suggested  with  the  use  of  spectral-holographic  four-wave  mixing  methods.  This  system  uses  a 
uiultipie-quantum  well  device  as  its  iionlliieaf  elenieiit.  At  this  point  It  is  not  clear  whether  the 
aberration,  quantum  1/f,  hysteresis,  and  device  non-uniformity  caused  noises  can  all  be  held 


sufficiently  low  at  the  MQW  device  in  order  to  make  this  suggestion  useful.  Refer  to  paper  #13  in 
the  list  of  publications  during  the  grant. 

6 .  The  new  concept  of  negative  quantum  entropy  has  been  successfully  applied  by  the  PI  to 
the  quantum  1/f  effect  for  the  first  time  in  1996,  and  presented  at  the  January  1997  WOFE 
workshop  in  Tenerife  ajid  at  the  May  1997  Orlando  Frequency  Control  Symposium.  The  concept 
had  been  developed  earlier  in  1996  and  explained  why  the  entropy  increases  in  measurements 
although  quantum  mechanical  evolution  is  unitary  and  precludes  entropy  increase  in  pure  states. 
The  Pi's  application  of  this  new  quantum  information-theoretical  concept  to  the  quantum  1/f  effect 
resolves  basic  questions  related  to  the  entropy  increase  generated  by  the  quantum  1/f  effect.  Refer 
to  papers  #^4  and  15  in  the  list  of  publications  during  the  grant. 

7 .  Finally,  the  application  of  the  quantum  i/f  theory  to  high-Q  quartz  resonators  and 
frequency  standards  mentioned  above  in  points  2  and  12  above  was  extended  for  the  first  time  to 
the  general  incoherent  case  of  lower-Q  and  surface-acoustic-wave  quartz  resonators.  This 
development,  important  for  practical  technical  applications  and  for  the  optimization  of  quartz 
devices  and  systems  was  presented  at  the  May  1997  Annual  Frequency  Control  Symposium  in 
Orlando.  Refer  to  paper  #15  in  the  list  of  publications  during  the  grant. 

8 .  Other  progress,  so  far  unpublished,  included  work  on:  1)  Quantum  1/f  formulas  for 
FET  and  HFET  Quantum  1/f  noise  optimization;  2)  Development  of  quantum  1/f  formulas  for  the 
three  independent  noise  sources  in  BJTs  and  HBTs;  3)  Methods  for  BIT  and  HBT  base  resistance 
evaluation  from  1/fnoise  measurements;  4)  Quantum  1/f  optimization  formulas  for  FET  and  HFET 
materials;  5)  Quantum  1/f  material  optimization  formulas  for  j  unction  type  devices:  j  unctions,  BJTs 
and  HBTs;  6)  Quantum  1/f  optimization  of  amplifiers;  7)  Quantum  1/f  optimization  of  oscillators. 

The  present  state  concerning  l/f  noise  in  high-tech  applications  is  still  characterized  by  a 
serious  discrepancy.  A  new  understanding  of  the  fundamental  quantum  nature  of  low-frequency 
(1/f)  noise  has  been  achieved,  but  in  practice,  most  high-technology  devices  still  operate  far  above 
the  quantum  limit  defined  by  the  new  theory.  Actually,  the  quantum  limit  can  be  achieved,  and 
even  this  limit  could  be  lowered  in  practice  by  improving  the  design  on  the  basis  of  the  quantum  1/f 
theory, 

The  trouble  is  localized  both  at  the  level  of  the  materials  used  and  at  the  level  of  device 
design.  Other  detrimental  aspects  are  the  presence  of  large  leakage  currents  in  heterostructure  field 
effect  transistors  (HFET,  studied  at  the  Wright  Patterson  AFB  by  F.  Schuermeyer)  or  junction 
devices  (studied  at  Hanscom  Rome  AFB  by  J.M.  Mooney  and  N.F.  Yannoni),  and  the  lack  of  a 
satisfactory  quantum  theory  of  mesoscopic  devices  which  starts  from  a  quantum  admittance  theory 
of  highly  transmissive  devices.  The  former  causes  very  detrimental  signal  power  drain  as  well  as 
large  i/f  noise,  while  the  latter  causes  device  builders  to  overlook  many  possibilities  offered  by 
quantum  interference  and  tunneling  effects  in  submicron  devices. 

The  quantum  1/f  effect  is  a  fundamental  fluctuation  of  ail  currents  of  charged  particles  j,  of 

all  scattering  cross  sections  o^,  recombination  cross  sections  o^,  tunneling  rates  or  transition  rates 
of  any  kind  F,  in  short  of  all  physical  cross  sections  o  and  process  rates  F,  given  by  the  universal 
formula  S(f)  =  2aA/fN  (conventional  quantum  1/f  equation  [1-13])  for  small  devices,  and  = 

2a/jtfN  (coherent  quantum  1/f  equation  [13-16])  for  large  devices.  These  two  forms  can  be 
combined  into  a  single  general  formula,  as  we  show  below.  Here  S(f)  is  the  spectral  density  of 

fractional  fluctuations  in  current,  6j/j,  in  scattering  or  recombination  cross  section  bo/o,  or  in  any 

other  process  rate  6F/F.  a=e2/hc=l/137  is  Sommerf eld's  fine  structure  constant,  a  magic  number 
of  our  world  depending  only  on  Planck's  constant  h,  the  charge  of  the  electron  e  and  the  speed  of 

light  in  vacuum  c.  A=2(Av/c)-/3jt  is  essentially  the  square  of  the  vector  velocity  change  Av  of  the 
scattered  particles  in  the  scattering  process  whose  fluctuations  we  are  considering,  in  units  of  c. 

Finally,  N  is  the  number  of  particles  used  to  define  the  notion  of  current  j,  of  cross  section  a  or  of 
process  rate  F. 
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II.  INCOHERENCE  OF  THE  QUANTUM  1/F  EFFECT  IN  LOW-Q  BAW  AND 

SAW  QUARTZ  RESONATORS 


1.  Introduction 


The  application  of  the  quantum  1/f  theory  to  bulk  acoustic  wave  (BAW)  and  surface 
acoustic  wave  (SAW)  quartz  resonators  has  so  far  been  limited  to  the  case  of  very  high  Q 
resonators  in  which  the  phonons  are  coherent  throughout  the  resonator  volume.  Here  the  quantum 
1/f  theory  is  applied  to  the  general  case  of  an  arbitrary  coherence  length  of  the  phonons,  which  may 
be  large  or  small  compared  with  the  size  of  the  quartz  resonator.  This  allows  to  extend  the  theory 
for  the  first  time  to  low-Q  resonators  in  which  the  phonons  are  localized  in  a  part  of  the  resonator 
volume.  The  theory  is  also  extended  to  include  defect  scattering  along  with  the  phonon  scattering 
case  exclusively  published  so  far. 

Phonon  scattering  in  the  resonator  is  known  to  limit  the  short-  and  medium-term  frequency 
stability  in  all  quartz  resonators  [  17J.  Phonon  scattering  can  occur  on  other  phonons  (particularly 
at  higher  temperatures)  or  on  crystal  defects  (favored  by  default  at  low  temperatures).  In  both 
cases  this  process  is  shown  to  yield  a  1/f  spectrum  of  resonator  frequency  fluctuations  through  the 
conventional  Quantum  1/f  nifect.  As  was  first  shown  on  this  basis  with  the  help  of  a  simple 
harmonic  oscillator  model  1 18],  BAW  and  SAW  quartz  resonators  ought  to  have  a  Q-^  dependence 
of  their  FM  power  spectrum.  1  his  has  been  experimentally  verified  by  Gagnepain  and  Uebersfeid 
for  BAW  resonators  [19]  when  they  first  noticed  their  1/Q4  4  |aw,  and  by  Parker  for  the  SAW 
case[20].  Although  the  quantum  1/f  effect  provided  the  historical  basis  for  the  derivation  of  the  Q-4 
law  [18]  as  being  caused  by  fluctuations  in  the  dissipation  rate  of  the  quartz  resonator,  the  exact 
mechanism  through  which  the  quantum  1/f  effect  modulates  the  dissipation  rate  remained  unknown 
from  1978  to  1991. 

Finally,  the  bridge  directly  connecting  l/f  noise  in  frequency  standards  to  the  quantum  1/f 

effect  was  discovered  [21  [-[23].  Here  is  how  it  works.  The  rate  F  of  photon-interactions  which 
remove  phonons  from  the  main  quartz  resonator  mode  is  modulated  by  the  quantum  1/f  effect, 
therefore  exhibiting  observable  quantum  fluctuations,  while  its  expectation  value  remains  constant. 
Indeed,  whenever  a  phonon  is  removed  from  the  main  resonator  mode,  the  time-derivative  of  the 

polarization  vector  of  the  quartz  crystal  dP/dt=P,  is  suddenly  jolted,  suffering  a  step-like 

modification  as  a  function  of  time.  From  Maxwell's  equations  we  know,  however,  that  P  is  added 
to  the  current  J  and  that  such  a  modification  of  the  current  causes  radiation.  Solving  Maxwell's 

equations  we  find  that  as  a  result  of  the  phonon  removal  there  is  a  constant  energy  of 

(l/43tEo)4(AP)-/3c3  radiated  away  per  unit  frequency,  i.e.  per  Hertz  at  any  frequency  f.  Dividing 

this  result  by  the  energy  of  a  photon  hf,  we  find  that  there  is  thus  a  probability  of  2a(AP)2/3jT:fe2c2 
for  the  emission  (radiation)  of  a  bremsstrahlung  photon  of  frequency  f.  Here  a  = 

(l/47t8o)(23te-/hc)  =  1/137  is  Sommerfeld's  fine  structure  constant,  a  dimensionless  universal 
constant  constructed  from  Planck's  constant,  the  charge  e  of  the  electron,  and  the  speed  of  light  c, 
SI  units  are  used  here,  while  Gaussian  units  were  used  in  [23]. 

Since  there  is  a  probability  of  2ci( AP)-/3jife-c-  «  1  for  the  emission  of  a  photon  of 
Irequency  f,  the  quartz  crystal  suifers  a  reaction,  or  a  recoil,  in  its  quantum  state,  causing  the 

phonon-emission  rate  F  to  perform  quantum  oscillations  with  frequency  f  and  with  two-sided 
spectral  density  S'  of  fractional  fluctuations  given  by  the  same  expression,  S'6r/r(f)  = 
2a(AP)-/3jtfe2c-.  This  is  the  quantum  1/f  effect.  The  one-sided  spectrum  is  thus 


7 


(1) 


S6r/r(f)  =  4a(  AP)2/3jtfe2c2. 


This  means  that  any  radiation  caused  or  implied  by  a  quantum  transition  from  one  state  to  another 
comes  with  a  price:  it  reacts  back  on  the  system,  causing  the  rate  of  that  transition  to  be  modulated 
by  exhibiting  observable  macroscopic  quantum  fluctuations  with  a  spectral  density  of  fractional  rate 
fluctuations  identical  to  the  photon  emission  probability  accompanying  the  transition  considered. 
No  knowledge  of  quantum  mechanics  is  therei^ore  needed  in  order  to  apply  the  quantum  1/f  effect, 
provided  one  manages  to  divide  the  energy  radiated  by  the  energy  hf  of  one  photon  of  frequency  f, 
i.e.,  provided  one  accepts  just  the  reality  of  Planck's  constant  h  and  Planck's  relation  between 
photon  energy  and  frequency.  Knowledge  of  electrodynamics  is  needed,  however,  in  order  to 
calculate  the  energy  radiated  in  a  transition. 

The  reader  interested  in  a  basic  understanding  of  the  quantum  1/f  effect  will  find  a  most 
accessible  description  at  the  end  of  p.  8  and  beginning  of  p.  9  in  [24].  That  description  considers 
scattering  of  electrons  as  an  example  of  transition  which  emits  radiation  and  suffers  a  quantum  1/f 
modulation  of  its  rate,  rather  than  considering  scattering  of  phonons,  which,  as  we  believe,  is  most 
important  in  quartz  resonators.  All  that  is  involved  in  that  derivation  is  the  notion  of  DeBroglie 
wave  associated  to  a  particle,  or  the  notion  of  wave  function.  Pre-quantum  mechanics  notions  are 
thus  sufficient  for  a  basic  understanding  of  the  recoil,  or  energy-loss  mechanism,  of  the  quantum 
1/f  effect  (Ql/fE).  For  a  practical  application  of  the  Ql/fE,  however,  classical  physics  is  sufficient. 


2.  Spatial  incoherence  of  phonon  loss  rate  fluctuations  in  low  and  intermediate  Q 
resonators. 


The  treatment  [21]-[23]  we  have  provided  so  far,  assumes  that  the  photons  are  spread  over 
the  whole  crystal,  and  that  therefore  there  is  coherence  of  the  quantum  1/f  phonon  loss  rate 
fluctuations  throughout  the  resonator  volume.  In  the  limit  of  very  high  Q  resonators,  and  of  high 

resonator  frequency,  the  mean  free  path  and  the  coherence  length  e  of  the  phonons  exceed  the  size 
of  the  resonator  and  our  assumption  is  Justified.  However,  in  the  low  Q  and  low  frequency  case 
the  coherence  length  is  small  compared  with  the  dimensions  of  the  resonator  crystal.  The  resonator 

volume  is  then  composed  of  many  incoherent  regions  of  volume  which  fluctuate  independently. 
This  is  applicable  in  particular  to  SAW  resonators,  as  suggested  empirically  by  Parker  et  al.  [25], 

[26].  We  are  here  borrowing  their  notation  of  the  size  of  the  coherent  volume  elements  by  e. 
Considering  the  v=y/e^  independently  fluctuating  regions  similar,  we  replace  Eq.  (1)  by 


S6r/r(f)=  y<(6ri/r)2>f  =  v-l<(6ri/ri)2>f  =4a(APi)2/3n:fve2c2 


(2) 


Here  we  assumed  that  r=vri  and  that  <(6ri/ri)->f  =  vSsr/rCO  is  independent  of  i.  With  v=V/£ 
we  finally  obtain 


Ssr/rCO  =4a83(APi)“/3n;fVe2c2  (3) 

in  the  incoherent  domain  i.e.,  for  sufficiently  large  V  and  small  Q.  The  corresponding  fluctuations 
in  the  frequency  cd  of  the  quartz  resonator  are  derived  from  the  equations 

(o2  =  a)o2-2r2,  (o6o)  =  -2r6r;  S6a)/(o(f)  =  (l/4Q4)S6v/Y(f),  (4) 

in  which  coq  would  be  the  natural  frequency  of  the  unloaded  quartz  resonator  mode  in  the  absence 
of  the  dissipation  T.  Therefore  the  spectral  density  of  fractional  frequency  fluctuations  Ssm/wIO 

8 


will  display  the  same  1/V  dependence  in  the  incoherent  regime,  down  to  a  volume  Below 

this  volume  we  expect  a  proportionality  of  S6a)/(o(f)  with  the  volume.  Consequently,  Sscj/ooCO 
will  first  increase  proportional  to  1/V  in  the  incoherent  region,  and  will  decrease  proportional  to  V 

when  V  is  lowered  to  values  below  e.  This  line  of  thought,  however,  neglects  the  fact  that 
resonators  operate  usually  in  the  lowest  vibration  mode,  and  that  therefore  a  decrease  in  volume 

will  correspond  to  a  certain  increase  in  the  frequency  Wo.  The  latter  is  also  connected  with  the 

speed  of  the  sound  Cs  and  the  acoustic  attenuation  length  I=l/aac=2csCy(0o*E- 

The  phonon  mean  free  path  is  about  40  A  for  bulk  wave  phonons  in  quatz  at  room 
temperature  and  500  A  at  liquid  nitrogen  temperatures.  This  approximates  the  phonon  coherence 

length  8  very  well.  For  SAW  phonons  the  corresponding  coherence  length  values  may  be  4  times 
lower  due  to  the  smaller  velocity  of  the  surface  wave  and  due  to  its  stronger  scattering. 

In  a  SAW  quartz  resonator  the  wave  is  localized  within  about  two  coherence  lengths  s  from 

the  surface.  Therefore,  in  the  incoherent  regime  V  =  28A  is  a  good  approximation. 
Consequently,  we  expect  an  increase  of  S6«)/a,(f)  proportional  with  1/A  when  the  resonant  area  A 
is  decreased,  down  to  very  small  areas  of  the  order  of  8-. 

The  incoherence  encountered  here  is  similar  to  the  incoherence  introduced  spatially  by  the 
veiy  small  coherence  length  of  electrons  in  semiconductors,  of  the  order  of  3oA.  The  quantum  1/f 
ettect  and  therefore  the  measured  I/f  noise,  has  no  spatial  correlations  in  semiconductors  down  to 
this  very  small  length  scale.  No  other  1/f  noise  "models"  can  explain  this  experimental  fact.  The 
quanturn  1/f  contributions  from  each  electron  are  independent,  which  causes  a  factor  i/Ne  to  be 
present  in  all  quantum  1/f  effect  formulas.  Here  Ne  is  the  number  of  electrons  present  in  the 
sample,  which  had  the  transition  considered  as  their  last  interaction. 

In  the  following  section  we  consider  the  quantum  1/f  effect  in  a  resonator  volume  of  size 

equal  or  smaller  than  8  in  all  directions,  the  result  being  applicable  whether  or  not  it  represents  the 
whole  resonator  mode. 


III.  SPATIAL  COHERENCE  OF  THE  Ql/f  PHONON  LOSS  RATE 
FLUCTUATIONS  IN  HIGH-Q  RESONATORS  INCLUDING  DEFECT 

SCATTERING. 

Let  the  resonator  volume  V  be  smaller  than  the  phonon  coherence  length  e  in  all  directions. 

Then  v=l  and  we  can  use  Eq.  (1).  In  Eqs.  (l)-(3),  (AP)-  is  the  square  of  the  dipole  moment  rate 
change  associated  with  the  process  causing  the  removal  of  a  phonon  from  the  main  oscillator 
mode:  scattering  of  a  main  resonator  mode  phonon  on  a  themal  phonon  of  higher  frequency 

<to>«kT/h.  After  this,  we  will  include  also  the  case  of  defect  scattering  below.  To  calculate 

(AP)-,  we  write  the  energy  W  of  the  interacting  mode  <(0>  in  the  form 

W  =  nfi<(o>  =  2(Nm/2)(dx/dt)- 

=(Nm/e2)(e  dx/dt)2=(m/Ne2)82(P)2;  (5) 

The  factor  two  includes  the  potential  energy  contribution.  Here  m  is  the  reduced  mass  of 
the  elementary  oscillating  dipoles,  e  their  charge,  g  a  polarization  constant  of  the  order  of  the  unity, 

and  N  their  number  in  the  resonator.  Applying  a  variation  An=l  we  get 

An/n  =  2IAPI/IPI,  or  AP=P/2n.  (6) 
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Solving  Eq.  (5)  for  P  and  substituting  into  (6),  we  obtain 

lAPI  =  (Nli<0)>/n)l/2(e/2g).  (7) 

Substituting  AP  into  Eq.  (1),  we  get 

r‘“Sr(f)  =  Nali<CL)>/3n:tmc2fg2  =  A/f  .  (8) 


This  result  is  applicable  to  the  fluctuations  in  the  loss  rate  F  of  the  resonator  volume. 

From  Eq.  (4)  the  corresponding  resonance  frequency  fluctuations  of  the  quartz  is  given  for 

V^£  by 


co“^S(o(f)  =  (l/4Q^)(A/0  =  Nafi<co>/12njtmc“fg“Q4;  (9) 

where  Q  =to/2r  is  the  quality  factor  of  the  single-mode  resonator  considered,  and  <a)>  is  not  the 
circular  frequency  of  the  main  resonator  mode,  but  rather  the  practically  constant  frequency  of  the 

average  interacting  (thermal)  phonon.  Indeed,  there  are  an  average  n(o=kT/fia)  phonons  present  in 
any  mode  of  frequency  (o.  For  the  case  of  quartz  resonators  we  have  used  the  interacting  thermal 
mode  of  average  frequency  <(o>  to  calculate  the  quantum  1/f  effect.  The  corresponding  AP  in  the 

main  resonator  mode  of  frequency  tOo  has  to  be  also  included,  but  is  negligible  because  of  the  very 
large  number  n  of  phonons  present  in  the  main  resonator  mode  and  entering  in  the  denominator  of 
Eqs.  (7)-(9). 

Considering  also  Eq.  (3),  Eq.  (9)  can  be  written  in  general  with  NaVN/V  in  the  form 

S(0  =  P'V/fQ4,forV<E,  (10) 

and 

S(f)  =  p'  £2/fVQ4,  for  V^8 ,  (11) 

where,  with  an  intermediary  value  <co>=10^/s,  with  n=kT/h<(0>,  T=3(X)K  and  kT=4  ‘lO'^l  J 
P'  =(N/V)ah<a»/12nn:g2mc2 

=1022(l/137)(10-27l08)2/i2kTjcl0-27  9  1020  (12) 

For  the  case  of  defect  scattering,  a  two-phonon  process  takes  place.  A  phonon  from  the 
main  resonator  mode  scatters  on  a  defect  and  a  phonon  of  comparable  frequency  emerges  into 

another  mode  with  much  smaller  phonon  occupation  number  n(o=kT/ho).  In  this  case  we  have  to 

replace <0)>  by  to  and  n<(o>  with  Uej,  which  gives  a  p-value  which  is  (<(0>/a))2  smaller,  i.e.  lO^- 

10^  times  smaller.  In  general,  therefore,  writing  r=r'  +  F",  we  obtain  for  the  combined  phonon 
and  defect  scattering  case,  in  general, 

P=P’  [F  2  +  (<a»/co)2F"2]/F2.  (13) 

Although  the  defect  scattering  term  is  small  at  room  temperature,  it  may  become  dominant  at  low 
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temperatures,  when  the  phonon  scattering  rate  F'  becomes  much  smaller  than  the  defect  scattering 
rate  F" . 

The  form  of  Eqs.  (10)-(13)  shows  that  the  level  of  1/f  frequency  noise  depends  not  only  as 
as  previously  proposed  [18],  but  also  on  the  oscillation  frequency  or  the  volume  of  the  active 
region.  This  theory  qualitatively  fits  the  data  of  Gagnepain  who  varied  the  Q-factor  with 
temperature  in  the  same  quartz  resonator  (but  not  frequency  or  volume),  the  data  of  Walls  who 
considered  several  quartz  reonators  which  differ  in  volume  and  frequency,  and  the  data  of  Parker 
for  SAW  resonators,  with  their  relatively  low  Q  values. 

The  theory  also  provides  the  basis  for  predicting  from  first  principles,  without  adjustable 
parameters,  how  to  improve  the  i/f  level  of  resonators,  beyond  just  improving  the  Q-factor,  which 
has  been  known  for  many  years,  and  which  has  been  related  [18]  to  fluctuations  in  the  dissipation. 
Since  the  1/f  noise  level  depends  on  the  active  volume,  in  the  coherent  regime  one  should  use  the 
lowest  overtone  and  smallest  diameter  consistent  with  other  circuit  parameters.  In  the  incoherent 
(low  Q)  case,  the  opposite  should  be  considered. 


IV.  QUANTUM  1/F  NOISE  IN  GaN  SAMPLES 

Due  to  the  importance  of  gallium  nitride  (GaN)  for  both  the  solid-state  microelectronics  and 
the  opto-electronics  of  the  future,  a  calculation  of  the  fundamental  1/f  noise  in  this  material  is  of 
primary  importance  and  is  long  overdue.  We  start  with  (1)  a  simple  evaluation  of  the  conventional 
quantum  1/f  effect  in  impurity  or  defect  scattering,  continue  with  (2)  the  calculation  of  the  same 
effect  in  acoustic  phonon  scattering  and  in  (3)  the  particularly  important  optical  phonon  scattering. 
On  this  basis  we  calculate  then  (4)  the  resulting  conventional  quantum  1/f  effect.  This  calculation 
will  be  done  both  for  n-type  and  p-type  devices  approximately,  by  neglecting  corrections 
introduced  by  the  energy  distribution  of  the  electrons  and  the  author's  cross-correlations  formula. 
The  errors  introduced  by  this  approximation  are  estimated  to  be  in  the  10-20%  interval.  We  then 
calculate(5)  the  coherent  quantum  1/f  effect,  and  provide  on  this  basis  (6)  a  final  result  for  GaN 
samples  of  any  size  and  any  nature. 

1.  Impurity  or  defect  scattering.  As  briefly  mentioned  at  the  end  of  Sec.  I  and  as  derived  below  in 

Sec.  V,  the  conventional  quantum  1/f  effect  in  a  scattering  cross  section  o  or  recombination  rate  F 
is  given  by 

S6a/o(f)  =  Ssr/KD  =  2aA/fN  (14) 

Here  S(D  is  the  spectral  density  of  fractional  fluctuations  in  current,  6j/j,  in  scattering  or 
recombination  cross  section  6a/a,  or  in  any  other  process  rate,  6F/F.  As  mentioned  above, 

a=e2/hc=l/137  is  Sommerfeld's  fine  structure  constant,  a  basic  number  of  our  world  depending 
only  on  Planck's  constant  h,  on  the  charge  of  the  electron  e  and  on  the  speed  of  light  in  vacuum  c. 
A=2(Av/c)2/33t  is  essentially  the  square  of  the  vector  velocity  change  Av  of  the  scattered  particles 
in  the  scattering,  recombination  or  tunneling  process  whose  fluctuations  we  are  considering,  in 
units  of  c.  Finally,  N  is  the  number  of  particles  used  to  define  the  notion  of  current  j,  of  cross 
section  o  or  of  process  rate  F. 

For  impurity  or  defect  -caused  scattering  we  obtain  therefore  the  quantum  1/f  coefficient 
2aA  =  (4a/3jt)(Av/c)-  =  (4a/33t)(Ap/meftx:)“  ,  (15) 

which  is  evaluated  assuming  a  thermal  energy  3kT/2  for  the  motion  of  the  current  carriers.  Note 
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that4a/3n  is  3.1’  10’^.  Thus, 

2aA  =  (4a/3oT)<(Av/c)2>  =  (4a/3jt)(6kT/niefrc2)  =  1.2- 10-9  (T/300K)(mo/mefr)  (16) 
and  we  obtain  for  electrons  with  effective  mass  mefiF  0.2mo  in  n-type  GaN 


2aA(e)  =6.]o-y  (17) 

and  for  holes  with  effective  mass  mefp=  0.8mo  in  p-type  GaN 

2aA(h)  =  1.5- 10-9.  (13) 

For  comparison,  we  note  that  for  n-type  GaAs  with  me(p=0.068mo  we  had  2aAe  =  1.8*  10-^  and 
for  p-type  GaAs  with  metr=0.5mo  we  had  2aAh  =  2.4- 10-9.  Here  mo  is  the  free  electron  mass. 


2-  Normal  phonon  scattering.  In  this  case  we  put  in  Eq.  (15)  lApI  =  AE/s,  where  <AE>=:kT  is  the 
energy  change  of  the  acoustical  phonon  in  the  scattering  process  and  s  is  the  speed  of  sound.  We 
obtain 

2aA  =  (4a/3re)(kT/mefrcs)2  =  3- 1 0-«  (T/300K)2(mo/meff)2.  ( 19) 

Therefore,  for  electrons  with  effective  mass  meip=  0.2mo  in  n-type  GaN 

2aA(e)  =  7.5- 10-7  (20) 

and  for  holes  with  effective  mass  meff=  0.8mo  in  p-type  GaN 

2aA(h)  =  4.7- 10-8.  (2i) 

3.  Optical  phonon  scattering.  In  this  case,  the  fractional  quantum  1/f  fluctuation  bT/Tof  the 

electron  scattering  rate  T  on  a  polar  optical  phonon  of  momentum  Rq  and  energy  R(0,^0)o  is 
described  by  the  quantum  1/f  coefficient 


2aA  =  (4a/3jt)<(Rq/mefrc)2>  =  (4a/3jt)2<(Rk/mefrc)2> 

=  (4a/3jt)(6kT/meiTc2)  =  1 .2- 10-9  (T/300K)(mo/meff) ,  (22) 

the  same  as  for  impurity  or  defect  scattering,  and  we  obtain  as  before  for  electrons  with  effective 
mass  mefr=  0.2mo  in  n-type  GaN 

2aA(e)  =  6- 10-9  (23) 

and  for  holes  with  effective  mass  mefF=  0.8mo  in  p-type  GaN 

2aA(h)  =  1.5*  10-9.  (24) 

For  comparison,  we  note  that  for  n-type  GaAs  with  mefr=0.068mo  we  had  2aA(e)  =  1.8*  10-8  and 
for  p-type  GaAs  with  me(r=0.5mo  we  had  2aA(i')  =  2.4*  10-9. 
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4.  Resulting  conventional  quantum  1/f  coefficient  and  spectral  density.  Both  for  electrons  and 
holes,  the  resulting  quantum  1/f  coefficient  of  fluctuations  in  the  mobility  ^  is  given  by  the  relation 

2aA(®>^)  =  2i(q,/fXi)2  2aAi(®’h)^  (25) 

where  2aAi  is  the  quantum  1/f  coefficient  calculated  for  the  scattering  process  #i  which  contributes 
to  limiting  the  observed  mobility  according  to  the  relation 

1/fA  =  I]i(l/pi)-  (26) 

Finally,  if  both  electrons  and  holes  contribute  to  the  conductivity  X,=e(p|Xh  +  npe).  the  spectral 
density  of  conventional  quantum  1/f  fractional  fluctuations  in  k  is 

S6X/?.(f)  =  2(e2a/  a2)[(pjih)2A(h)/Nh  +  (nHe)2A(e)/Ne] .  (27) 

This  allows  for  the  calculation  of  quantum  1/f  noise  in  any  small  sample.  The  transition  to  the 
coherent  quantum  1/f  effect  applicable  for  large  samples  is  discussed  at  the  end  of  the  next  Sec.  V. 

Our  general  conclusion  is  encouraging,  and  can  only  help  to  brighten  the  future  for  GaN 
even  more.  Specifically,  as  we  have  seen  above,  the  conclusion  is  that  small  GaN  samples  and 
devices  offer  lower  1/f  noise  levels  than  comparable  GaAs  samples  or  devices,  by  a  factor  of  3  to 
10,  mainly  due  to  the  larger  effective  masses  of  the  current  carriers. 


V.  NEGATIVE  ENTROPY  IN  THE  QUANTUM  1/f  EFFECT 
1.  Schematic  derivation  of  the  conventional  quantum  1/f  effect 

The  author's  recent  application  of  the  new  Quantum  Information  Theory  Approach  (QIT, 
developed  1996)  to  Infra  Quantum  Physics  (IQP)  explains  for  the  first  time  the  apparent  lack  of 
unitarity  caused  by  the  entropy  increase  in  the  Quantum  1/f  Effect  (Ql/fE,  developed  by  us  1974). 
Indeed,  he  now  offers  a  more  rigorous  proof  of  the  conventional  quantum  1/f  effect  in  this  report, 
which  shows  that  actually  there  is  no  resultant  entropy  increase  and  therefore  unitarity  is  not 
violated.  His  new  proof  involves  the  concept  of  von  Neumann  Quantum  Entropy,  including  the 
negative  conditional  entropy  concept  for  quantum  entangled  states.  The  Ql/fE  was  applied  to 
many  high-tech  systems,  in  particular  to  ultra-small  electronic  devices.  The  present  report  explains 
how  the  additional  entropy  implied  by  the  Ql/fE  arises  in  spite  of  the  entropy-conserving  evolution 
of  the  system. 

The  conventional  quantum  1/f  effect  is  a  fundamental  fluctuation  of  all  currents  of  charged 
particles),  of  all  scattering  cross  sections  recombination  cross  sections  a^.,  tunneling  rates  or 

transition  rates  of  any  kind  T,  in  short  of  all  physical  cross  sections  o  and  process  rates  T,  given 
by  the  universal  formula 

S(f)  -  2aA/lN  (28) 

We  present  here  only  a  schematic  derivation.  Let's  simplify  our  world  and  assume  only 
one  electromagnetic  mode  of  the  universe  would  be  present,  with  frequency  co  and  wave  vector  k. 
Consider  the  field  mode  in  its  ground  state  and  a  pair  of  2  incoming  identical  charged  particles  with 
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the  same  well-defined  wave  vector  being  scattered  by  some  potential.  Both  the  initial  and  the  final 
state  represent  a  pure  state.  The  initial  state  is  l++->o,  where  we  reserved  the  first  two  arguments 
of  the  ket  for  the  two  electrons  and  the  last  (-)  for  the  field.  Due  to  the  interaction  with  the  field, 
the  final  state  will  be 

|f>  =  ([++->  +  y\+-+>/V2  +  y\-++>/V2)/iV\+y'^),  (29) 

where  y  is  the  emission  amplitude  of  a  bremsstrahlung  photon,  i.e.,  for  the  excitation  of  the  field 
oscillator  from  its  ground  state  (-)  to  its  first  excited  state  (-I-).  The  first  two  arguments  label  the 
state  of  the  charged  particles,  indicating  the  presence  of  an  energy  loss  with  (-)  and  the  persistence 
in  the  same  energy  state  with  (+).  The  third  argument  of  the  kets  always  labels  the  field  oscillator, 
as  mentioned. 

The  corresponding  density  operator  of  the  pure  state  obtained  is  p  =  IfxfI.  Its  quantum 
von  Neumann  entropy  S/k  =  a  =  -Tr(p  log  p)  is  zero. 

2.  Paradoxical  entropy  increase  in  1/f  noise  and  quantum  information  theory 

Ignoring  the  field  oscillator,  i.e.,  taking  the  trace  of  over  the  field  oscillator  label,  we  obtain  a 
classically  correlated  system,  a  mixture  of  2  pure  states,  described  by  the  density  operator 

[I++X++1  +  yH\+->  +  l-+>)(<4-l  +  <-+l)/2]/(l-hY2).  (30) 

The  second  term  gives  the  QI/fE  at  f=(o/2n:,  as  we  show  below.  The  corresponding  entropy  is 

log(l+Y^)  -  (Y^logY^)/(l+Y^)  >  0  for O^^^l.  (31) 

The  entropy  of  the  system  thus  appears  to  have  increased  although  according  to  quantum 
mechanics  it  can  not  increase.  Indeed,  according  to  quantum  mechanics,  time  evolution  of  a  state 

occurs  through  a  unitary  transformation.  The  latter,  however,  is  known  to  leave  a  =  -Tr(p  log  p) 
invariant. 

The  quantum  information  theory  (QIT)  solves  this  paradox.  When  two  systems  A  and  B 
with  quantum  von  Neumann  entropy  o(A)  =  -Tr^(p^^logp^)  and  o(B)  =  -Trjj(pjjlogpjj)  form  a 
composite  system  AB  of  entropy  o(AB)  =  -Tr^(p^logp^Q)  we  can  prove  that  we  have  to  write 

a(AB)  =  a(A)  +  a(BIA)  =  o(B)  +  o(AIB),  (32) 

in  perfect  analogy  with  classical  entropies  or  information.  We  have  introduced  o(AIB)  as  the  von 
Neumann  entropy  of  A  conditional  on  B,  or  the  entropy  of  A  when  we  know  B: 

a(AIB)  =  ■Tr^[p^logp^uj] .  (33) 

Here  we  have  introduced  the  conditional  density  matrix  p^|g  =  Paj3(1a®  Pb)’^’  quantum  analog 
similar  to  the  classical  conditional  probability,  where  @  is  the  tensor  product  in  the  joint  Hilbert 

space  and  marginal  density  matrix  obtained  by  taking  a  partial  trace  over  the 

variables  associated  with  A. 

3.  Negative  entropy  and  solution  of  the  paradox 

The  conditional  entropy  is  usually  negative  in  quantum  entangled  states.  Finally,  we 
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introduce  the  quantum  mutual  entropy  which  represents  the  shared  entropy,  corresponding  to  the 
mutual  information  between  A  and  B: 

o(A:B)  =  -Tr[p,^logp^  3]  s  a(A)  +  o(B)  -  o(AB),  (34) 

where 

Pa:B=Pab(Pa®Pb)‘-  (35) 

Taking  all  logarithms  in  the  base  of  2,  the  entropies  will  be  expressed  in  bites.  Considering  our 

paradoxical  case  discussed  above  for  simplicity  with  y^=1.  we  obtain  the  following  entropy 
diagram  of  our  quantum  mechanically  entangled  triolet  of  three  systems  comprising  the  charged 
particles  A  and  B.  as  well  as  the  field  oscillator  C: 


Fig.  1. 

Here  the  quantum  conditional  entropies  are  underlined  for  emphasis,  and  the  mutual  quantum 
entropy  has  double  underlining.  Dotted  underlining  marks  the  mutual  entropy  of  two  of  the 
subsystems,  which  is,  however,  not  the  mutual  entropy  of  the  whole  system.  We  see  that  A,  B 
and  C  have  each  1  bit  of  entropy,  as  we  expect,  since  each  can  be  in  a  +  or  -  state.  In  addition, 

®ABC  =  0  expected  for  a  pure  state.  The  negative  conditional  entropies  indicate  that  this  state  is  a 
pure  quantum  state  which  can  not  be  obtained  in  classical  physics. 

If  we  ignore  C  by  tracing  over  it,  we  obtain  the  classically  correlated  system  AB  with  its 
positive  entropy  of  1  and  the  negative  quantum  entropy  state  of  the  ignored  field  oscillator, 
conditional  on  AB: 


In  conclusion,  the  negative  conditional  entropy  of  the  ignored  field  oscillator 
compensates  the  positive  entropy  of  the  system  of  two  particles.  Unitarity  and  entropy 
conservation  are  both  satisfied.  This  is  the  solution  of  the  paradox.  It  explains  why  our  earlier 
explanation  of  the  Ql/fE  in  terms  of  a  two-particle  wave  function  was  correct,  in  spite  of  the 
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apparent  lack  of  unitarity. 

4.  Derivation  of  the  conventional  anantum  1/f  effect 

In  the  rest  of  the  present  section  we  apply  the  scheme  developed  here  to  the  real  world  with 
all  field  oscillators  present  in  l++->o  and  with  N  incident  particles  in  l++->o. 

We  start  with  the  expression  of  the  Heisenberg  representation  state  if>  of  N  identical 

bosons  of  mass  M  emerging  at  an  angle  6  from  some  scattering  process  with  undetermined 

bremsstrahlung  energy  losses  reflected  in  their  one-particle  waves  tpi(^i) 

lf>  =  (N!)-l/2  Hi /d3|i(pi(|,)tp+(ii)  I0>  =  Hi /d3|,(pi(§i)  |p>>,  (36) 

where  is  the  field  operator  creating  a  boson  with  position  vector  and  I0>  is  the  vacuum 

state,  while  lfo>  is  the  vacuum  field  state  with  N  bosons  of  position  vectors  |i  with  i  =  All 

products  and  sums  in  this  section  run  from  1  to  N,  unless  otherwise  stated. 

To  calculate  the  particle  density  autocorrelation  function  in  the  outgoing  scattered  wave,  we 
need  the  expectation  value  of  the  operator 

0(xi,X2)  =  '»l>+(xi>tlt+(x2H>(x2>il)(xi),  (37) 

known  as  the  operator  of  the  pair  correlation.  Using  the  commutation  properties  of  the  boson  field 
operators,  we  first  calculate  the  matrix  element 

N!<So|Olf‘b.= 

2  l)^('np“*2)^(?n~Xl)^(^m— X2)2](i  j)n'ij6(T)j— §i).  (38) 

Here  the  prime  excludes  p=v  and  m=n  in  the  summations  and  excludes  i=m,  i=n,  j=(ji  and  j=v  in 

the  product.  The  summation  2(i  j)  runs  over  all  permutations  of  the  remaining  N-2  values  of  i  and 
j.  On  this  basis  we  now  calculate  the  complete  matrix  element 

<flOlf>=  [1/N(N-1)] 

Tp  (^p)T\'  (^v)Tm(im)Tn(in)^('nv~*l)^(llp'*2)^(^n"Xl)^(?m‘*2)'  (39) 

=  [1/N(N-1)1  S'^v^'mn  Cpp*(X2)<Pv*(xi)(Pm(xi)cpn(X2) 

The  one-particle  states  are  spherical  waves  emerging  from  the  scattering  center  located  at 

x=0: 


cp(x)  =  (C/x)e>i^’‘  [1  +  a+k,i].  (40) 

Here  C  is  an  amplitude  factor,  K  the  boson  wave  vector  magnitude,  y  =  b(k,l)  the  bremsstrahlung 

amplitude  for  photons  of  wave  vector  k  and  polarization  1,  while  a'*'k,i  is  the  corresponding  photon 
creation  operator,  allowing  the  photon  state  to  be  created  from  the  vacuum  if  Eq.  (40)  is  inserted 

into  Eq.  (36).  The  momentum  magnitude  loss  q=Mck/hKsMf/hK  is  necessary  for  energy 
conservation  in  the  Bremsstrahlung  process.  Substituting  Eq,  (10)  into  Eq.  (9),  we  obtain 

<flOlf>  =  IC/x|4{N(N-1)  +  2(N-l)2k,ilb(k,l)|2f  l+cosq(xi-X2)]},  (41) 

where  we  neglected  a  small  term  of  higher  order  in  b(k,l).  To  perform  the  angular  part  of  the 
summation  in  Eq.  (41),  we  calculate  the  current  expectation  value  of  the  state  in  Eq.  (40),  and 
compare  it  to  the  well  known  cross  section  without  and  with  bremsstrahlung 

j  =  (hK/Mx2)[l  +2kllb(k,l)|2]=jjl  +/aAdf/f],  (42) 
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where  the  quantum  fluctuations  have  disappeared,  where  a=e^/hc  is  the  fine  structure  constant, 

aA=(2a/37i)(Av/c)-  is  the  fractional  bremsstrahlung  rate  coefficient,  also  known  in  QED  as  the 
infrared  exponent,  and  the  1/f  dependence  of  the  bremsstrahlung  part  displays  the  well-known 
infrared  catastrophe,  i.e.,  the  emission  of  a  logarithmically  divergent  number  of  photons  in  the  low 

frequency  limit.  Here  Av  is  the  velocity  change  h(K-Ko)/M  of  the  scattered  boson,  and  f=ck/2jt 
the  photon  frequency.  Eq.  (41)  gives 

<flOlf>  =  IC/x|4{N(N-1)  +  2(N-l)(JaAdf/f)[l+cosq(xi-X2)]},  (43) 

which  is  the  pair  correlation  function,  or  density  autocorrelation  function  along  the  scattered  beam 
with  df/f=dq/q.  The  spatial  distribution  fluctuations  along  the  scattered  beam  will  also  be  observed 
as  fluctuations  in  time  at  the  detector,  at  any  frequency  f.  According  to  the  Wiener-Khintchine 

theorem,  we  obtain  the  spectral  density  of  fractional  scattered  particle  density  p,  (or  current  j,  or 

cross  section  o)  fluctuations  in  frequency  f  or  wave  number  q  by  dividing  the  coefficient  of  the 
cosine  by  the  constant  term  N(N-l): 

p-2Sp(f)  =  j-2Sj(f)  =  cr2So(f)  =  2aA/fN  [or  =2aA/f(N-l)  for  fermions],  (44) 

where  N  is  the  number  of  particles  or  current  carriers  used  to  define  the  current  j  whose 
fluctuations  we  are  studying.  Quantum  1/f  noise  is  thus  a  fundamental  1/N  effect. 

I'he  exact  value  of  the  exponent  of  f  in  Eq.  (44)  can  be  determined  by  including  the 

contributions  from  all  real  and  virtual  multiphoton  processes  of  any  order,  and  turns  out  to  be  aA- 

1,  rather  than  -1,  which  is  important  only  philosophically,  since  aA«l.  The  spectral  integral  is 
thus  convergent  at  f=0. 

For  fermions  we  repeat  the  calculation  replacing  in  the  derivation  of  Eq.  (38)  the 
commutators  of  field  operators  by  anticommutators,  which  finally  yields  in  the  same  way 


p-“Sp(f)  =  j-2Sj(f)  =  Gr2s^f)  =  2aA/f(N-l),  (45) 

which  causes  no  difficulties,  since  N>2  for  particle  correlations  to  be  defined,  and  which  is 
practically  the  same  as  Eq.  (44),  since  usually  N»l.  Eqs.  (44)  and  (45)  suggest  a  new  notion  of 
physical  cross  sections  and  process  rates  which  contain  1/f  noise,  and  express  a  fundamental  law 
of  physics,  important  in  most  high-tech  applications. 

S,  Derivation  of  the  coherent  quantum  1/f  effect 

The  present  derivation  is  based  on  the  well-known  propagator  Gs(x'-x)  derived 
relativisti(^lly  [27-31]  in  a  new  QED  picture  required  by  the  infinite  range  of  the  Coulomb 
potential.  1  he  corresponding  nonrelativistic  form  was  provided  by  Zhang  and  Handel  [32]: 


-i<0onril)s<x')ipj(x)l<l>o>  =  6ss'  Gs(x'-x) 

=  (i/V)2{expi[p(r-r')-p2(t-t’)/2m]/ll}np^s 

P 

x{-ip(r-r’)/h+i(m2c2-i-p2)l/2(t-t')(c/h)}a/3i.  (46) 

Here  a  is  again  Sommerfeld's  fine  structure  constant,  Up  s  the  number  of  electrons  in  the  state  of 

momentum  p  and  spin  s,  m  the  rest  mass  of  the  fermions,  6ss'  the  Kronecker  symbol,  c  the  speed 
of  light,  x=(r,t)  any  space-time  point  and  V  the  volume  of  a  normalization  box.  T  is  the  time- 
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ordering  operator  which  orders  the  operators  in  the  order  of  decreasing  times  from  left  to  right  and 
multiplies  the  result  by  (-1)^,  where  P  is  the  parity  of  the  permutation  required  to  achieve  this 
order.  For  equal  times,  T  normal-orders  the  operators,  i.e.,  for  t=t'  the  left-hand  side  of  Eq.  (46) 

is  i<<bohl>J(x)i|)s'(x')l‘J)o>.  The  state  <I>o  of  the  N  electrons  is  described  by  a  Slater  determinant  of 
single-particle  orbitals.  Eq.  (46)  is  exact  at  large  t-t',  or  low  frequency. 

To  calculate  the  current  autocorrelation  function  we  need  the  density  correlation  function, 
which  is  also  known  as  the  two-particle  correlation  function.  The  two-particle  correlation  function 
is  defined  by 

<d>orri|)f(x)%(x>ipKx')il)s<x')10o>  =  <d>oltl>J(x)ilJs(x)l<I>o><^ol^?(x')il’s<x')l<l>o> 

-  <d)orrips<x')ipit(x)l<I>(P><OoTtl)s(x)'ipKx')IOo>.  (47) 

The  first  term  can  be  expressed  in  terms  of  the  particle  density  of  spin  s,  n/2  =  N/2V  = 

<^olipf(x)'il)s(x)l^tP>-,  while  the  second  term  can  be  expressed  in  terms  of  the  Green  function  (22) 
in  the  form 

Ass'(x-x')  S<d>oltl)j(x)ipj'(x')'»ps<x')tt)s(x)l4)o>  =  (n/2)-  +6ss'  Gs(x'-x)Gs(x-x').  (48) 

The  "relative"  autocorrelation  function  A(x-x')  describing  the  normalized  pair  correlation 
independent  of  spin  is  obtained  by  dividing  by  n-  and  summing  over  s  and  s' 

A(x-x')  =  1  -  (l/n2)2Gs(x-x')Gs(x’-x) 

S 

=  1  -  (1/N2)2  5;{expi[(p-p')(r-r’)-(Ep-EpO(t-t')]/ft}np,snp',s 

s  pp' 

x{p(r-r')/h-(m-c--i-p-)^/-(t-t')(c/h)}“/^ 

x{p’  (r'  -r)/fi-(m2c2-hp'  -)i/2(t'-t)(c/R)}«/K.  (49) 

Here  we  have  used  Eq.  (46).  We  now  consider  a  beam  of  charged  fermions,  e.g.,  electrons, 

represented  in  momentum  space  by  a  sphere  of  radius  pp,  centered  on  the  momentum  po  which  is 
the  average  momentum  of  the  fermions.  I'he  energy  and  momentum  differences  between  terms  of 
different  p  are  large,  leading  to  rapid  oscillations  in  space  and  time  which  contain  only  high- 
frequency  quantum  fluctuations.  The  low-frequency  and  low-wavenumber  part  A]  of  this  relative 
density  autocorrelation  function  is  given  by  the  terms  with  p  =  p' 

Ai(x-x')  =  1  -  (1/N2)2  S^p.s 

s  p 

xlp(r-r')/R-(m-e-+p-)T2(t-t')(c/li)l2«/’i  (50) 

«  1  -  (l/N)lpo(r-r’)/R-mc2T/R|2«/’t  for  pp«lpo3-mc2T/zl.  (51) 

Here  we  have  used  the  mean  value  theorem,  considering  the  2a/jt  power  as  a  slowly  varying 

function  of  p  and  neglecting  po  in  the  coefficient  of  x  s  t-t',  with  z  s  Ir-r'I.  The  correlations 
propagate  along  the  beam  with  a  group  velocity  given  by  the  average  velocity  po/m  of  the  particles 

in  the  beam,  and  with  the  phase  velocity  of  c-N.  Denoting  6  s  It  -po(r-r')/mc-l,  and  using  a  well 
known  identity  for  the  last  step  [32],  we  obtain  from  Eq.  (51)  the  form 

Ai(x-x')  =  1  -  (1/N)lmc20/R)|2a/Jt  =  i  -(1 .25/N)l0P«/’i  =  1  -(1.25/N)e(2a/3t)ln0 
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«  l-(1.25/N)[]+(2a/jt)ln01  =  1  -  (2.5/N)  +  (1.25/N)[l-(2a/jt)ln0] 


(52) 


«(N-2.5)/N  +(1.25/N)e-(2«/Jt)lne^  (iyN)(N-2.5+(2.5a/:tcosg)  fcos(o0da)/(ol-2«/^}. 


This  indicates  a  or  1/f  spectrum  and  a  l/(N-2.5)  dependence  of  the  spectrum  of  fractional 

fluctuations  in  density  n  and  current  j.  The  total  error  corresponding  to  the  two  linear 
approximations  of  exponentials  in  Eq.  (52)  is  less  than  1%,  provided  lln0l<2O,  or  (250,000)’* 
hours  <  101  <  250,000  hours.  Here  0=6/ls  and  to  is  the  circular  Fourier  frequency  in  rad/s.  We 

have  used  [Is  mc-ZK)]-"^’*  «  1 .25;  this  accounts  also  for  the  presence  of  the  number  2.5  instead  of 
the  more  normal  number  2  in  the  final  form.  The  form  we  have  chosen  here  is  more  convenient  for 
applications.  The  equivalent  normal  form  would  have  been 

00 

Ai(x-x')  « (l/N){N-2+(2a/jtcosa)  J[mc“/lia)  ]-“^^  costoO  dto/w},  (53) 

in  which  the  error  caused  by  the  two  linear  approximations  of  exponentials  would  have  been  of  the 
order  of  20%,  and  in  which  the  fractional  power  would  also  have  been  neglected  in  the  integrand 

for  all  purposes  except  for  the  theoretical  question  of  the  integrability  of  the  l/(0  spectrum.  The 
fractional  autocorrelation  of  current  fluctuations  6j  is  obtained  by  multiplying  Eqs.  (52-53)  on  both 
sides  with  (epo/m)-,  and  dividing  by  (enp(ym)-  which  is  the  square  of  the  average  current  density 
j,  instead  of  just  dividing  by  n-.  So  it  is  the  same  as  the  fractional  autocorrelation  for  quantum 
density  fluctuations.  The  last  form  of  Eq.  (52)  for  the  universal  coherent  quantum- 
electrodynamical  chaos  process  in  electric  currents  finally  becomes 

S6j/j(to)  «  [2.5a/jtto(N-2.5)][co]2«/’t  «  2.5a/n;ajN  =  0.0058/coN.  (54) 

The  spectral  density  resulting  from  Eq.  (53)  coincides  with  the  result  2a/jtfN,  derived 
directly  earlier  [14]  from  the  coherent  state  of  the  electromagnetic  field  of  a  physical  charged 
particle.  The  connection  with  the  conventional  quantum  1/f  effect  was  discussed  in  Sec.  II.  1. 

Being  observed  in  condensed  matter  in  the  presence  of  a  constant  applied  field,  these 
fundamental  quantum  current  fluctuations  are  usually  interpreted  as  mobility  fluctuations.  Most  of 
the  conventional  quantum  1/f  fluctuations  in  physical  cross  sections  and  process  rates  are  also 
mobility  fluctuations,  but  some  are  also  in  the  recombination  speed  or  tunneling  rate. 

6.  Discussion  and  conclusions 


Consider  now  the  physical  basis  of  the  connection  with  the  conventional  quantum  1/f 
effect.  The  energy  of  the  collective  (drift)  motion  of  the  current  carriers  has  a  magnetical 
component,  which  we  consider  coherent,  because  the  field  contributions  of  the  individual  current 
carriers  are  addi  ti  ve 

Em  =  /(B2/8:r)d3x  =  [nevS/c]2ln(R/r).  (55) 

Here  n  is  the  concentration  of  current  carriers,  v  their  drift  velocity,  c  the  speed  of  light,  r  the 
radius  and  S  the  area  of  the  cross  section  of  the  cylindrical  conductor  or  semiconductor  sample 
considered.  1  here  is,  however,  also  a  kinetic  energy  contribution,  which  is  incoherent,  because 
the  kinetic  energies  of  the  drift  motion  of  the  current  carriers  must  be  added  to  get  the  total  kinetic 
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energy,  rather  than  their  momenta. 


Ej|(^  =  Emv2/2  =  nSmv2/2  =  Ejjj/s.  (56) 

The  "coherece  ratio" 

s  =  Ej^/E]j  =  2ne2s/mc2ln  (R/r) »  2e2N7mc2,  (57) 

has  been  introduced  here  for  convenience,  where  N'  =  nS  is  the  number  of  carriers  per  unit  length 
of  the  sample  and  the  natural  logarithm  ln(R/r)  has  been  approximated  by  one  in  the  last  form.  We 
expect  the  observed  spectral  density  of  the  mobility  fluctuations  to  be  given  by  a  relation  of  the 
form 


(l/fx2)S^(f)  =  [l/(l+s)][2aA/fN]  +  [s/(l+s)][2a/jtfN]  (58) 

which  can  be  interpreted  as  an  expression  of  the  effective  Hooge  constant  if  the  number  N  of 
carriers  in  the  (homogeneous)  sample  is  brought  to  the  numerator  of  the  left  hand  side 

ttH  =  [l/(l+s)]aHconv  +  [s/(  l+s)]aHcoh-  (59) 

In  this  equation  aA=2a(Av/c)2/3n:  is  the  usual  nonrelativistic  expression  of  the  infrared 
exponent,  present  in  the  familiar  form  of  the  conventional  quantum  1/f  effect.  This  equation  is 
limited  to  quantum  1/f  mobility  (or  diffusion)  fluctuations,  and  does  not  include  the  quantum  1/f 
noise  in  the  surface  and  bulk  recombination  cross  sections,  in  the  surface  and  bulk  trapping 
centers,  in  tunneling  and  injection  processes,  in  emission  or  in  transitions  between  two  solids. 

Note  that  the  coherence  ratio  s  introduced  here  equals  the  unity  for  the  critical  value  N'  = 
N"  =  2*10l2/cm.,  e.g.  for  a  cross  section  S  =  2*10"'4  cm^  of  the  sample  when  n  =  10^^.  For 
small  samples  with  N'«N"  only  the  first  term  survives,  while  for  N'>N”  the  second  term  in  Eq. 
(59)  is  dominant. 

The  conventional  and  coherent  Ql/fE  was  successfully  applied  FETs,  HFETs,  BJTs,  and 
HBTs,  as  well  as  to  quartz  resonators  and  amplifiers.  In  conclusion: 

—  The  quantum  1/f  effect  appears  to  generate  entropy,  but  this  is  only  because  we  do  not 
observe  the  negative  conditional  entropy  contribution  of  the  emitted  soft  photons. 

—  The  quantum  1/f  effect  contains  Planck's  constant  in  the  denominator  and  diverges  in  the 
classical  limit. 

“  The  quantum  1/f  effect  represents  a  new  aspect  of  quantum  physics. 

—  The  quantum  1/f  effect  is  a  form  of  quantum  chaos,  arising  only  from  spontaneous 
bremsstrahlung. 

—  The  quantum  /f  effect  limits  most  high-technology  applications  and  devices.  Its 
knowledge  allows  us  to  optimize  them. 


VI.  CHARACTERISTIC  FUNCTIONAL  OF  PHYSICAL  THERMAL  NOISE 
-  WHICH  INCLUDES  EQUILIBRIUM  QUANTUM  1/f  NOISE  - 

1.  Introduction 


In  a  previous  paper  [33]  the  characteristic  functional  of  quantum  1/f  noise  was  derived. 
This  result  was  applicable  to  quantum  1/f  noise  in  any  cross  section  or  process  rate,  and  in  currents 
or  voltages.  In  another  paper  [34]  the  application  of  quantum  l/f  noise  to  thermal  equilibrium 
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noise  was  performed,  and  the  quantum  1/f  Nyquist  theorem  was  formulated.  Deviations  from  the 
Gaussian  form  of  thermal  noise  were  expressed  [34], [35]  in  terms  of  the  characteristic  function. 

The  Quantum  1/f  Effect  (Ql/fE)  is  a  fundamental  fluctuation  phenomenon  present  as 

coherent  Ql/f  Noise  2a/jtfN  in  electrical  currents  of  any  nature  and  present  as  conventional  Ql/f 


Noise  2aA/fN  in  the  physical  cross  sections  and  process  rates  (PCS&PR).  The  PCS&PR  are 
defined  as  the  scattered  current,  or  the  outgoing  current  resulting  from  a  process,  including  its 
quantum  fluctuations  acquired  due  to  the  interaction  with  the  electromagnetic  field,  all  referred  to 
the  unit  of  incoming  flux,  which  is  one  particle  per  second  per  cm-.  The  usual  notion  of  cross 
section  or  process  rate  is  a  quantum-mechanical  expectation  value  which  does  not  contain  the 
quantum  fluctuations.  In  the  present  report  the  nature  of  the  quantum  1/f  fluctuations  introduced 
by  us  in  1975  is  investigated  and  described  as  a  form  of  quantum  chaos. 

Both  Coherent  and  Conventional  Quantum  1/f  Noise  are  peculiar  quantum  effects.  They 
contain  Planck's  constant  only  in  the  denominator,  as  part  of  the  fine  structure  constant  to  which 
both  are  proportional.  This  implies  that  they  increase  without  upper  bound  when  we  approach  the 
classical  limit  by  letting  h  go  to  zero.  One  is  left  with  the  impression  that  the  classical  limit  of  the 
quantum  1/f  effects  does  not  exist. 

At  this  point  it  is  useful  to  remember  that  the  quantum  1/f  approach  originated  from  the 
classical  homogeneous  and  isotropic  turbulence  theory  in  an  unbounded  plasma  of  electrons  and 
holes  [36],  [37],  and  was  introduced  as  the  result  of  prolonged  quantization  efforts  applied  to  this 
classical  theory.  The  classical  theory  corresponded  to  infinite  noise  power  in  any  finite  frequency 
interval.  Therefore,  if  we  accept  quantum  1/f  noise  as  the  quantum  manifestation  of  classical 
turbulence,  its  infinite  classical  limit  no  longer  confronts  us  with  conceptual  difficulties. 
Moreover,  the  physical  nature  of  quantum  1/f  noise  becomes  connected  directly  to  the  nature  of 
classical  turbulence  which  is  in  practice  caused  by  the  instability  of  laminar  flow  when  a  certain 
critical  value  of  the  Reynolds  number  is  surpassed. 

Being  governed  by  the  Navier-Stokes  equations,  classical  turbulence,  naturally  obtained 
from  the  instability  of  laminar  flows,  is  in  fact  deterministic,  in  spite  of  its  stochastic  appearance. 
Therefore  we  use  the  word  "chaos"  to  describe  its  true  nature.  The  transition  from  laminar  flow  to 
turbulence  is  one  of  the  seven  known  routes  to  chaos,  the  oldest  known  route  in  fact.  The  chaotic 
nature  of  turbulence  characterizes  turbulence  as  a  deterministic  process  if  finite  dimensionality 
dependent  on  the  bandwidth  of  the  turbulent  frequencies  observed.  In  the  same  way,  quantum  1/f 
noise,  being  governed  by  the  Schrodinger  Equation  for  the  chaotic  wave  fields,  is  in  fact 
deterministic;  it  represents  the  oldest  known  form  of  quantum  chaos.  The  finite  dimension  of 
quantum  1/f  noise  was  calculated  [38]  on  the  basis  of  a  1 -dimensional  model  of  the  quantum  1/f 
theory. 


In  conclusion,  quantum  chaos  is  defined  here  as  the  quantum  manifestation  of  classical 
chaos.  The  latter,  in  turn,  is  defined  as  the  quasistationary  state  eventually  reached  by  a  system 
after  instabilities  of  the  laminar  state  have  triggered  the  transition  into  another  state  characterized  by 
time-dependent  values  of  the  dynamical  parameters,  with  a  probability  distribution  law  which  is 
stable  in  time.  Our  notion  of  quantum  chaos  is  visualized  by  the  deterministic  evolution  of  the 
seemingly  stochastic  Schrodinger  field. 

As  we  know,  however,  the  wave  function  is  not  directly  observable,  but  rather  is  a  means 
of  calculating  the  expectation  values  of  observable  quantities.  Even  if  the  wave  function  is  well 
defined,  due  to  its  probabilistic  interpretation,  we  obtain  random,  or  stochastic  fluctuations  of  the 
local  physical  parameters  of  the  system,  superposed  on  the  deterministic-chaotic  background.  We 
will  consider  this  stochastic  component  irrelevant  for  the  following  reason:  the  stochastic  element 
enters  only  because  we  ask  for  properties  which  the  system  does  not  have,  i.e.,  because  we  asked 
a  meaningless  question.  In  quantum  mechanics  such  meaningless  questions  can  be  made 
artificially  meaningful  by  imposing  the  corresponding  properties  on  the  system  from  outside  in  the 
measurement  process.  This  does  not  mean  the  system  had  this  property  all  along,  it  had  instead 
other,  cornpleiiientary,  paranieters  which  were  better  defined  before,  but  lost  their  definition  in  the 
measuring  process.  The  physical  meaning  of  the  superposed  stochastic  component  is  that  it 
represents  the  shot  noise,  which  is  well  known,  given  by  the  spectral  density  of  2el,  and  which  is 
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not  the  object  of  our  studies.  It  is  worth  mentioning,  however,  that  it  is  a  organical  component  of 
the  noise  process  which  has  a  deterministic  (chaotic),  as  well  as  a  superposed  stochastic 
component. 

The  (31/fE  is  caused  by  the  infinite  range  of  the  electromagnetic  interaction,  i.e.,  of  the 
Coulomb  field.  Any  charged  particle  is  thus  always  in  interaction  with  all  other  arbitrarily  distant 
charges,  and  this  causes  the  particle  to  have  an  energy,  which  is  not  well  defined,  i.e.,  a  mass,  and 
therefore  a  relativistic  mass  shell,  which  is  not  sharp.  Mathematically,  this  energy  uncertainty  is 
expressed  by  the  fact  that  the  charged  particle  consists  of  the  bare  particle  plus  its  electromagnetic 
field  which  in  the  rest  system  is  just  the  Coulomb  field.  The  latter  is  characterized  by  a  well- 
defined  electric  field  and  vector  potential.  A  Fourier  transformation  then  tells  us  that  the 

electromagnetic  field  oscillators  for  any  wave  vector  k  and  polarization  X  have  a  well  defined 
phase,  and  are  therefore  in  coherent  states.  The  latter,  in  turn,  are  states  of  indefinite  energy,  first 
introduced  by  Schrodinger  in  the  early  days  of  quantum  mechanics,  and  later  used  by  Glauber  to 
describe  coherent  states  of  the  electromagnetic  field. 

How  can  we  express  mathematically  the  stochastic  fluctuations  and  their  deterministic- 
chaotic  background?  Heisenberg's  equations  of  motion  for  the  density  matrix  of  the  system  of 
charged  particles  whose  quantum  1/f  current  fluctuations  are  under  consideration  and  of  the 
electromagnetic  field  modes  coupled  to  them  provide  the  desired  description.  The  solution  should 
be  given  by  the  characteristic  functional  of  the  current  fluctuations.  The  characteristic  functional  of 
conventional  quantum  1/f  noise  j(t)  =  J(t)  -  <J>  affecting  the  current  I=<J>  has  been  shown  [6]  to  be 

H[x(8)]  =  exp[-aA(l2/N)/(8/eo)“^lx(e)|2d8/£] .  (60) 

Here  aA  is  the  infrared  exponent  aA=(4a/3jt)<(Av/c)2>,  which  is  known  to  enter  the  expression 
of  both  the  basic  quantum  1/f  noise  formula  and  the  characteristic  functional  of  quantum  1/f  noise 
in  two  ways:  as  a  coefficient  and  as  an  additional  exponent  of  the  frequency  f.  Sommerfeld's  fine 

structure  constant  a=e-/hc  =1/137  is  well  known,  and  so  is  the  Boltzmann  constant  ke-  Eq.  (60) 
gives  the  characteristic  functional  of  1/f  noise  for  a  sample  containing  N  scattered  current  carriers 

with  a  velocity  change  <(Av)“>  caused  the  scattering  process.  Eq.  (60)  also  approximates  the 
characteristic  functional  of  1/f  noise  from  a  solid-state  sample  containing  N  current  carriers  with  an 

average  velocity  change  <(Av)“>  in  the  scattering  processes  which  determine  their  mobility. 

The  present  report  performs  the  final  step  in  the  study  of  thermal  noise  with  infrared 
radiative  corrections,  by  calculating  the  characteristic  functional  of  physical  thermal  noise.  We  call 
it  "physical"  because  we  reler  to  the  actually  observed  thermal  equilibrium  (Johnson)  noise,  rather 
than  to  the  strictly  Gaussian  noise  expected  from  the  Nyquist  formula  without  infrared  radiative 
corrections.  The  small  quantum  1/f  contributions  come  in  as  time-dependent  infrared  radiative 
corrections,  required  by  the  interaction  of  the  current  carriers  with  the  electromagnetic  field,  or  by 
the  reaction  of  the  bremsstrahlung  back  on  the  current  which  has  produced  it.  Knowing  the 
characteristic  functional  of  a  process  is  the  best  we  can  do,  the  highest  level  and  most 
comprehensive  and  complete  description  possible  for  a  random  process. 

2«  Amplitudie  Distribution  and  Characteristic  Function 

In  terms  of  current  fluctuations  in  equilibrium,  to  derive  the  characteristic  function  of  the 

physical  thermal  noise  current  variable  |,  we  express  it  in  terms  of  the  unmodulated  theoretical 
Nyquist  noise  current  variable  x: 

I  =  x(G/Go)1^2  =  X  +  xy,  (61) 

where  y=6G/2Go  is  the  fractional  quantum  1/f  fluctuation  in  the  conductivity  G=Go-i-6G  of  the 
conductor  whose  thermal  equilibrium  current  fluctuations  we  are  considering.  Let  the  quantum  1/f 
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noise  variable  y  obey  a  Gaussian  amplitude  distribution  of  dispersion  02,  and  x  one  of  dispersion 
oi.  We  prove  now  that  in  terms  of  the  zero-order  Bessel  function  of  imaginary  argument,  the 

product  z=xy  will  have  an  amplitude  distribution  P(z)=(l/3taia2)Ko(z/oia2)  which  has  an 
elementary  characteristic  function 

X(v)  =  (1  +  (62) 

In  order  to  derive  the  amplitude  distribution  of  z  =  xy,  we  introduce  r|  =  In  z,  a  =  In  x  and 
P  =  In  y.  From  =  a  +  p  we  obtain  the  distribution 
P"(il)dTi  =  dr\fP'i(a)P'2(Ti-a)da 

=  (2dr|/2jT0i02)/exp[-(e2«)/20i^+a]exp[-e2(il-«)/202^+(Ti-a)]da 
00 

=(e'n  dTi/2jt0i02)  Jexp{-[e^  lo\02]  coshE}d8  =(e^  dTi/jtOi02)Ko(eWoi02),  (63) 

-00 

where  P'(oi)  and  P'2(P)  are  derived  from  Gaussians  of  dispersion  Oi  and  02  respectively,  e  =  2a  - 
IT  +  ln(a2/oi),  and  K  is  the  modified  zero-order  Bessel  function.  Returning  to  the  real  variable  z, 
we  obtain  from  P"('n)dT|  =  P(z)dz 

P(z)  =  (l/3iaia2)Ko(z/aia2).  (64) 

A  Fourier  transform  shows  that  the  corresponding  characteristic  function  is  elementary: 

Z(v)  =  (1  +  ai2a22v2)-l/2.  (55) 

This  is  just  what  we  wanted  to  prove. 

Our  next  step  is  the  calculation  of  the  characteristic  function  of  | .  Due  to  the  mutual 

dependence  of  x  and  z,  the  characteristic  function  of  |  differs  from  the  product  of  the  characteristic 
functions  of  x  and  z.  We  obtain  instead  the  amplitude  p.d.f. 

00 

I(i)=  fPi(x)P2(|/x-l)(l/lxl)dx 
-00 


00 

=  [exp(-l/2a2^)/2jtaia2]  Jexp[-x2/2ai2  -  ^'^12x^02^  +  |/xo2^]dx/lxl 

-00 

00 

=  [exp(-l/2a2^)/23raia2]^expl-x2/2ai2  -|2/2x2a22]cosh(|/xo22)dx/x,  (66) 

which  we  have  written  down  directly,  but  which  can  also  be  obtained  by  substituting  (e“  -  1)~  for 

6^“  in  Eq.  (63).  This  distribution  tends  to  a  Gaussian  in  the  limit  02-^0,  but  differs  slightly  from 
a  Gaussian  in  general.  The  skewness  is  zero,  but  the  kurtosis  exhibits  a  small  deviation  from  the 
Gaussian  amplitude  p.d.f.. 

By  performing  a  Fourier  transform  of  Eq.  (66)  we  find  again  that  the  characteristic  function 
of  the  total  physical  thermal  noise  is  elementary 

X(k)  =  (1  +  exp[-k2ai 2/2(1  +  k20i2022)].  (67) 

while  the  amplitude  distribution  itself  can  not  be  expressed  in  elementary  functions.  Here  0i2= 
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4kTGB  and  02^=(<xAJ2N)J{£lEo)*^^de/t],  with  B=F-f  representing  the  bandwidth  considered. 

In  the  remaining  part  of  this  section  we  derive  the  characteristic  functional  of  thermal  noise 
including  thermal  equilibrium  quantum  1/f  noise  in  the  conductance  G.  We  first  familiarize  the 
reader  with  characteristic  functionals  in  general. 

3.  Characteristic  Functional 


The  complete  characterization  of  a  random  process  i(t)  is  given  by  its  characteristic 
functional  [39] 

H[x(t)]=<ei(j.x)>.  (68) 

where  j(t)  is  in  our  case  the  current  (or  voltage,  or  fractional  conductance)  fluctuation  j(t)  =  J(t)  - 
<J>,  x(t)  is  the  argument  function  of  the  functional,  and  the  scalar  product  in  Hilbert  space  is 
defined  as 

T 

(j,x)  =  jj*(t)x(t)dt.  (69) 

Here  T  is  the  interval  of  observation  of  the  fluctuations,  i.e.,  the  length  of  the  sample  considered 
from  the  stationary  random  process  at  hand.  Particularly  simple,  the  characteristic  functional  of  a 
Gaussian  process  is  [39] 

HG[x(t)]  =  exp[-(l/2)<(j,x)2>].  (70) 

In  particular,  for  thermal  short-circuit  current  fluctuations  in  a  resistor  of  conductance  G 

described  by  the  simple  Nyquist  formula  with  oi^  =  4kTGB,  the  characteristic  functional  is 

F  T 

=  exp[-2kBTGBjlx(t)|2dt],  (71) 

where  Xf  is  the  Fourier-transformed  argument  function.  Here  the  two-sided  frequency  integral  is 
limited  to  the  allowed  bandwidth  B,  and  the  time  integral  to  the  interval  T  of  observation.  Since 

B— >00,  keT  must  be  replaced  by  hf/[e*’f^*^^-  1]  and  included  in  the  frequency  integral,  while  the 
time  integral  must  also  include  the  Fourier-transformed  Planck  kernel. 

The  characteristic  functional  of  quantum  1/f  noise  affecting  the  conductance  G  has  been 
shown  [6]  to  be 

H[x(8)]  =  exp[-aA(G2/N)j(e/8o)“^tx(e)|2d8/8] .  (72) 

Here  aA  is  the  infrared  exponent  aA=(2a/3n:)<(Av/c)->,  which  is  known  to  enter  the  expression 
of  both  the  basic  quantum  1/f  noise  formula  and  the  characteristic  functional  of  quantum  1/f  noise 
in  two  ways:  as  a  coefficient  and  as  an  additional  exponent  of  the  frequency  f.  Sommerfeld's  fine 

structure  constant  a=e-/hc  =1/137  is  well  known,  and  so  is  the  Boltzmann  constant  ks.  Eq.  (72) 
gives  the  characteristic  functional  of  1/f  noise  for  a  sample  containing  N  current  carriers  with  an 

average  velocity  change  <(A v  )2>  in  the  scattering  processes  which  determine  their  mobility.  Since 

we  are  interested  in  the  variable  y=6G/Go,  we  must  replace  G  by  1/4  in  Eq.  (72). 

For  each  value  of  the  frequency  we  now  use  Eqs.  (71)  and  (72),  repeating  the  steps  which 
led  to  the  derivation  of  Eq.  (65).  This  yields  the  characteristic  functional  of  the  term  z=xy  in  Eq. 
(61) 


H[x]  =  exp[-2GJkBTIx{l2df] 
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Z[v(f)]  =exp[-(l/2)/ln(l  +  Oi^a2^v^)d{\ 

=  exp{-(l/2)/ln[l  +  2kBTG(aAyfN)(f/fo)«A  v2(f)]df}.  (73) 

According  to  Eq.  (61)  we  now  have  to  add  the  Nyquist  term  x.  Repeating  for  every  frequency  f 
the  steps  which  led  us  to  Eq.  (67),  we  obtain  now  the  characteristic  functional  of  physical  thermal 

noise  ^  in  the  form 

F[k(f)]  =  exp{-(l/2)Jln[l  +  2kBTG(aA/fN)(f/fo)«Ak2(f)]df} 
exp{-2kBTG/k2(f)[l  +2kBTG(aA/fN)(f/fo)«Ak2(f)]-ldf}.  (74) 

Eq.  (74)  gives  the  characteristic  functional  of  physical  thermal  noise  for  a  sample  containing  N 
current  carriers  with  an  average  velocity  change  <(Av)2>  in  the  scattering  processes  which 
determine  their  mobility,  resulting  in  a  well  defined  infrared  exponent  aA. 


VII.  A  BRIDGE  BETWEEN  COHERENT  AND  CONVENTIONAL 

QUANTUM  1/f  NOISE 


In  this  section,  a  certain  mass  distribution  is  defined  which  allows  arbitrary  physical 
electronic  propagators  with  infrared  radiative  corrections  of  any  order  included  corresponding  to 
given  Hamiltonians  and  boundary  conditions  to  be  approximated  by  an  integral  over  the  product  of 
this  mass  distribution  with  the  corresponding  well-known  electronic  propagator  without  infrared 
radiative  corrections. 

Coherent  Quantum  1/f  Noise  Theory  [14]  was  developed  in  1982-83,  nine  years  after 
Conventional  Quantum  1/f  Noise  [1]  was  introduced.  Both  arc  forms  of  the  same  infrared 
divergence  phenomenon  and  are  described  by  simple  universal  formulas  involving  Sommerfeld's 

fine  structure  constant  a=e2/hc,  but  the  coherent  Ql/f  formula  2a/nf  is  particularly  simple. 
Conventional  Ql/f  Noise  is  always  present  in  any  physical  process  rate  or  cross  section,  but  in 
devices  which  are  sufficiently  large,  the  larger  coherent  Ql/f  effect  dominates  the  1/f  current  noise. 
Since  the  conventional  Ql/f  effect  is  associated  with  individual  carriers  and  the  coherent  Ql/f  noise 
with  their  collective  drift  motion,  a  physical  interpolation  formula  [15]  was  developed  in  1985, 
superposing  the  two  results  in  the  same  proportion  as  the  kinetic  and  magnetic  energies  of  the  drift 
motion  of  the  carriers. 

As  was  pointed  out  [13]  (See  also  the  end  of  Sec.  V  above),  the  coherent  Ql/f  effect  can 
also  be  derived  from  a  quantum-electrodynamic  propagator,  but  both  this  and  the  1983  derivation 
are  applicable  only  for  a  sufficiently  large  sample  or  device  in  which  the  energy  of  the  collective 
magnetic  field  of  the  drift  motion  far  exceeds  the  sum  of  the  additional  kinetic  energies  of  the 
carriers  due  to  their  drift.  For  smaller  samples  we  will  attempt  in  this  report  to  construct  a 
generalization  of  coherent  Ql/f  noise  which  is  based  on  the  physical  interpretation  of  a  charged 
particle  as  a  quantum  object  with  indefinite  energy  or  mass,  specifically,  with  a  fuzzy  mass  shell. 

To  achieve  this  objective,  we  represent  the  non-relativistic  form  [32]  of  the  new  quantum- 
electrodynamic  propagator  as  a  superposition  of  classical  propagators,  defined  by  an  unknown 


mass  distribution  p((x)  which  describes  the  fuzzy  mass  shell: 

■  2  *  2 
r??!lrv  -U  f  ^2  .  tixpa/ji 

iv-vi-r ti 


00 


Jdpp(|x).exp{  {y[v<r-r')-(c2+y)(t-t’)]}} .  (75) 
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1  ,  V- 

Let  u  =  j^[v.(r-r’)-(c^+^)(t-t')].  This  allows  us  to  simplify  the  above  equation,  in  the  form 

(76) 


00 


When  we  use  p'=  p-m,  the  equation  becomes 


00 


Jcip'p'(p')-e't^'“=(imu)«/’». 


(77) 


Because  p'(p')  is  different  from  zero  only  around  p'=0  or  p=m  ,  we  can  extend  the  domain  of 
integration. 


00 


J’dp'  p'  (p'  ).e’^^'  “=  (imu)®^’'. 


(78) 


-00 


Let  us  take  the  derivative  with  respect  to  u.  This  yields 


00 


a 


j*dp’p'(p’)^*^‘'“.ip’  =  ^ 


-00 


-a/jc 


(79) 


We  can  further  simplify  the  above  equation  with  the  notation  p*  .p'  (p')=  X(p*)  and  get 


00 


a 


*m 


,a/ji 


^dp'X(p')«‘^^'‘'=  -- 
-TO  (iu 


(iu)i-a/«  ■ 

We  can  determine  X(p')  by  taking  the  Fourier  transformation  of  the  right  hand  side. 


(80) 


a 


,a/jt 


00 


00 

X(p’)  =  fdu  — - =  — *5—  rdu 

-o)  (iu)‘-“'’^  2jt-i’-«'".i,  u‘-“'" 

_  am“^”  /  °  cos(p'u)+isin(p'u)  °°  cos(p'u)+isin(p'u)\ 

2ji2jl  a/ji  yl  a/z  J  a/z 

f  y.x  cos(-p'u')+isin(-p*u')  ”  cos(p'u)+isin(p■u)^ 
—  f^u’  -<^Qs(p’u’)+isin(p’u’)  ”  cos(p'u)+isin(p’u)x 
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(81) 
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(83) 


(for  fx'  >  0),  and 

r»  f‘TT 


f/t  .  n»-«/.)nctM) 


(for  fi'  <  0). 

Because  both  1 and  sin(a/2)  are  much  smaller  than  l-(-l)^'“^"^  and  cos(a/2),  we  can  just 
use 


X(u')  =(i  (  i)i-a/^)  ar(a/:t)cos(a/2)  m.g/^ 

2jr2jl-a/n 


(84) 


for  all  practical  purposes.  We  thus  conclude  that  the  mass  distribution  function  has  to  be 


P(1a)  = 


ar(a/jT)cos(a/2)  m“'" 


(85) 


This  is  a  remarkable  result.  It  allows  us  to  approximate  the  effect  of  infrared  radiative 
corrections  on  any  electronic  propagator  by  multiplying  it  by  p(p)  and  integrating  over  p  as  was 
done  with  the  free  particle  propagator  on  the  right-hand  side  of  our  first  equation  above.  The  result 
will  represent  an  approximation  of  the  physical  electron's  propagator  corresponding  to  the  problem 
at  hand,  i.e.,  an  approximation  of  the  physical  propagator  including  the  infrared  radiative 
coirections,  which  corresponds  to  the  given  potential  in  which  the  electron  has  to  move,  and  which 
satisfies  the  given  boundary  conditions. 


VII.  APPLICATION  OF  HIGH  e  PARAMAGNETIC  AND  FERROELECTRIC 
MATERIALS  TO  HIGH-TRANSCONDUCTANCE  HFETs  WITH 
LOW  NOISE  AND  LOW  GATE  CURRENT 

1.  Introduction 


The  main  problem  in  high-transconductance  FETs  is  the  presence  of  a  considerable  gate 
leakage  current  that  loads  the  input  and  reduces  the  power  amplification  coefficient.  The  problem 
has  become  particularly  serious  in  InGaAs  Heterostructure  FETs. 

Indeed,  to  achieve  a  high  transconductance,  the  gate  capacitance  must  be  high,  and  the  gate 
insulation  must  be  thin.  This  leads  to  tunneling  currents  through  the  gate  insulation. 

A  detailed  examination  of  the  quantum-mechanical  tunneling  problem  by  these  authors  in 
cooperation  with  F.  Schuermeyer  in  1993-1996  shows  that  the  only  way  to  reduce  the  gate 
tunneling  current  without  sacrificing  the  gate  capacitance  is  to  use  materials  with  high  dielectric 

constant  Er  as  insulators  under  the  gate. 

It  is  obvious  that  only  amorphous  ferroelectric  materials  with  negligible  remanent 
polarization  and  high  coercive  field  are  useful  if  we  want  to  avoid  the  hysteresis  effect,  and  are 
interested  only  in  suppression  of  the  leakage  currents  in  high  transconductance  devices.  For 
certain  applications,  hysteresis  in  the  response  of  the  device  is  desirable,  and  should  be  achievable 
by  slight  modifications  in  the  fabrication  process.  The  present  report  analyzes  the  possibility  of 

using  ferroelectric  and  paraelectric  high  8  gate  insulators,  including  also  noise  considerations,  in 
particular  quantum  1/f  noise. 

2.  General  Considerations 

An  increased  gate  capacitance  is  particularly  desirable  at  lower  frequencies,  because  it  leads 
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to  higher  transconductances  of  the  FET,  without  compromising  its  desirable  large  input  impedance. 
At  high  and  ultrahigh  frequencies,  however,  a  very  large  gate  capacitance  C  is  less  desirable, 
because  it  essentially  shorts  the  gate  to  the  channel  and  requires  large  capacitive  gate  currents  and 
signal  power  input  levels  at  the  achievable  fj  and  fmax  values.  The  natural  solution  that  allows  for 

a  smaller  frequency  dependence  of  the  input  admittance  jwC  is  a  ferroelectric  gate  insulation. 
Indeed,  at  low  frequencies,  it  could  increase  the  gate  capacitance  and  device  transconductance  by  a 
factor  of  10^,  while  at  higher  frequencies,  where  such  a  large  capacitance  could  be  harmful,  the 
dielectric  constant  decreases  to  intermediate  values,  allowing  for  operation  at  the  highest 
frequencies  with  a  reasonable  input  impedance.  This  is  the  main  advantage  of  a  ferroelectric  gate 
insulation. 

Unfortunately,  the  ferroelectric  introduces  its  own  resonances  and  their  characteristic 
frequencies,  as  well  as  dissipation.  Furthermore,  if  leakage  through  the  ferroelectric  is  still 
occurring,  it  leads  to  a  normalized  quantum  1/f  noise  coefficient  larger  by  a  factor  of  10^  than  what 
we  would  get  without  piezoelectric  coupling  of  the  current  carriers. 

A  possible  way  for  reducing  the  remaining  frequency  dependence  of  the  gate  capacitance  is 
the  use  of  two  or  more  ferroelectric  layers.  This  introduces  additional  resonances  and  resonances 
corresponding  to  the  space-charge  relaxation  at  the  interfaces.  However,  in  return,  each  of  these 
resonances  could  be  weighted  down  so  that  the  resulting  frequency  dependence  will  just  present 
small  nonuniformities  which  cause  negligible  spectral  distortion  in  amplifiers. 

3.  History 


The  effect  of  ferroelectric  polarization  on  insulated-gate  thin  film  transistor  parameters  was 
first  studied  by  Zuleeg  and  Wieder  [40]  as  a  device  for  non-volatile  memory  design.  ^They  applied 
a  1  nim  thick  film  of  ferroelectric  TGS  (tri-glicine-sulphate)  on  one  side  of  a  5,OOOA  CdS  film  as 
gate  insulation,  the  other  side  of  the  film  being  insulated  by  lOOOA  of  Si02  from  an  A1  gate.  Gold 
was  used  for  the  gate  applied  to  the  TGS,  the  source  and  the  drain  contacts.  Depending  on  the 
polarization  state  of  the  TGS,  controlled  by  the  gold  gate  voltage,  and  by  its  previous  voltage 
History,  One  gcis  various  sets  of  drain  characteristics.  In  the  nietal-ferroelectnc  semiconductor 


structure,  proposed  later  by  Wu  [41],  a  thin  film  of  ferroelectric  Bi4Ti30j2  was  replacing  the  oxide 
of  a  MOSFET.  This  was  suggested  as  a  new  ferroelectric  memory  devfce.  The  materials  PZT, 
i.e.,  PbZrxTii-xQj  [42]  and  LiNbOs  [43]  were  investigated  as  ferroelectric  gate  insulators,  as  well 
as  other  materials  [44j.  However,  all  these  materials  contain  oxygen,  and  will  cause  oxidation  of 
the  compound  semiconductors. 

un  the  other  hand,  the  nature  of  the  leakage  current  was  investigated  in  HFETs  by 
Schuermeyer  et.  al.,  [45]-[46].  They  focused  primarily  on  the  thermal  activation  of  carriers  from 
the  channel  over  the  ban  ier  separating  them  from  the  gate.  These  authors  found  that  the  density  of 
the  thermal  activation  current  reached  a  maximum  at  an  intermediary  region,  closer  to  the  drain  than 

to  the  source.  In  the  present  report,  we  illustrate  the  application  of  ferroelectric  or  high  e  materials 
to  the  case  of  the  InGaAs/AlGaAs/GaAs  and  InGaAs/InAlAs/InP  enhancement-mode  HFETs 
fabricated  and  studied  by  Schuermeyer  et.  al.,  [45]-[46]. 


4.  Approach 


We  will  consider  the  ferroelectric  material  BaMgF4  as  a  leading  candidate  for  coating  the 
barrier  layer,  applying  it  in  its  amorphous  form  at  450o  C  on  the  external  AlGaAs  barrier  layer, 
separating  the  channel  from  the  gate  of  the  HFET.  In  this  form  it  has  negligible  hysteresis  and  is 

characterized  by  a  dielectric  constant  £i=35.  The  deposition  temperature  is  not  too  high  for  the 
HFET.  The  AlG^s  barrier  thickness  d  can  be  reduced  to  15  nm  and  the  BaMgF4  layer  needs  to 
be  only  20  nm  thick,  f  ig.  l*and  2’show  the  approximate  device  structure  and  the  conduction  band 
configuration  for  an  n-channel  HFET. 

To  lurther  clarify  our  approach,  we  use  the  simple  FET  expression  for  the  intrinsic 
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Fig.  f:  HFET  structure. 


Fig.  2*:  HFET  Band  Structure 
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Here  Cq  is  the  gate  capacitance,  tc  the  effective  channel  transit  time,  e  the  effective  permittivity  of 
the  gate  capacitance,  L  the  gate  length,  W  its  width,  tins  the  effective  distance  between  gate  and 

channel,  equivalent  with  a  total  insulation  thickness,  p  the  effective  channel  mobility,  Vp  the 
source-drain  voltage,  and  Vm  the  maximal  drift  velocity  obtained  in  the  high  field  regime  of  the 
channel. 

From  Eqs.  (86)  and  (87)  we  see  that  by  increasing  tins  6y  a  factor  X  and  e  by  a  factor 
gm  will  increase  by  a  factor  of  X^.  At  the  same  time  the  leakage  current  caused  by  tunneling 
through  the  barrier  will  become  negligible  for  all  practical  purposes  for  X^2,  because  the  tunneling 

probability  T«  1  will  be  replaced  by  Indeed,  an  additional  decrease  of  the  gate  current  will  be 

caused  by  the  reduction  of  the  electric  field  present  throughout  the  barrier. 

The  unit  current  gain  frequency,  ft,  is  not  affected,  being  approximated  by 

ft  =  l/2jrtc.  (88) 

However,  the  maximum  frequency  of  oscillation,  f^ax*  is  given  by 
fmax  =  (fx/2)(gni/gout)^^^>  (89) 

where  gout  is  the  output  conductance  that  is  not  affected  by  Cq  or  by  X.  The  maximum  frequency 
of  oscillation  is  therefore  increased  by  provided  there  is  no  decrease  in  8  at  f=fmax-  For 
instance,  with  X=1.41,  increasing  tins  1-41  times  and  8  4  times  will  increase  gm  2.82  times,  fmax 
1.65  times  and  stop  the  gate  current.  Actually,  there  will  be  a  considerable  decrease  of  8  at  high 
frequencies,  and  therefore,  we  expect  fmax  to  increase  less  than  predicted  by  the  factor  X^^-. 

5.  BaMeFa  Technology 

The  ferroelectricity  of  BaMgF4  was  discovered  in  1969.  However,  this  material  was  not 
fabricated  as  a  thin  film  till  1989,  and  was  first  deposited  on  Si  surfaces  to  generate  ferroelectric 
memory  FBI's  [47].  It  has  been  found  that  BaMgF4  grown  on  GaAs  or  AlGaAs  (100)  surfaces  at 
450°C  is  amorphous  [48].  This  film  does  not  exhibit  hysteresis,  but  is  a  paraelectric  insulator  with 
a  relative  dielectric  constant  8r  of  35.  The  substrate  was  (100)  n-type  GaAs  at  450°C.  The  growth 
rate  was  0.4  pm/h. 

At  temperatures  of 500-600^^0  polycrystalline  ferroelectric  films  exhibiting  hysteresis  were 
grown,  with  a  (140)  orientation.  Their  coercive  field  was  about  200k V/cm  and  the  remanent 

polarization  Pr  about  1.3  pC/cm-.  For  monocrystalline  material  the  expected  value  of  Pr  is  4 
pC/cm^. 

In  order  to  obtain  high-gm  HFETs  with  truly  negligible  gate  currents,  we  therefore 
recommend  choosing  the  factor  X  defined  in  Sec.  4  to  be  v^2.  We  thus  obtain  a  transconductance 

increased  by  a  factor  of  2V2,  a  total  barrier  thickness  increased  by  a  factor  of  y/2,  and  therefore  a 
negligible  gate  current.  Practically,  this  is  achieved  by  properly  structuring  the  barrier  grown  on 
the  undoped  GaAs  channel.  First  a  5  nm  undoped  AlGaAs  spacer  layer  is  deposited  on  the 
channel,  then  10  nm  of  n-type  AlGaAs,  followed  by  20  nm  of  amorphous  BaMgF4,  and  by  the 
metal  of  the  gate  electrode,  e.g.,  Al. 

By  varying  the  growing  temperature  and  speed  in  the  interval  of 450-495OC  and  around  0.4 
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^im/h,  we  may  be  able  to  achieve  a  useful  further  increase  in  e  without  detrimental  hysteresis 
effects. 

6.  1/f  Noise 

The  1/f  noise  of  the  2-DEGrees-of-freedom-electrons  in  the  conducting  InGaAs  channel  is 
determined  by  the  conventional  quantum  1/f  formula  for  polar  optical  scattering  at  300K,  or  for  the 
mixture  of  p>olar  optical  phonon,  longitudinal  acoustical  phonon,  piezoelectric,  and  lattice  defect 

scattering  at  70K.  The  spectral  density  S  of  fractional  quantum  1/f  fluctuations  ST/T  of  the  electron 

scattering  rate  T  in  a  process  with  velocity  change  Av  of  the  electrons  is 

S6r/r(f)  =  2aA/fN  =  (4a/3jtfN)(Av/c)2.  (90) 

The  fractional  quantum  1/f  fluctuation  bT/Tof  the  electron  scattering  rate  T  on  a  polar  optical 
phonon  of  momentum  hq  and  energy  htOq^cOo  is  given  by  the  spectral  density 

Ssr/rCO  =  (4a/3jtfN)<(hq/meffc)2>  =  (4a/33tfN)2<(fik/meffc)2>.  (91) 

Here  we  have  estimated  the  average  quadratic  momentum  change  of  the  electrons,  resulting  from 
the  momentum  fiq  of  the  optical  phonon  to  be  2(fik:)2,  because  the  momentum  transfer  to  electrons 
of  momentum  fik  is  known  [46]  to  be  between  the  limits  hqniin=  fi{-k+[k2  +  2meffC0o/R]^^2].  and 

h{k+[k2  +  2meflt0o/H]^^2}..  With  <(hk/meiT)“>  =3kT/mefl'  and  mefp  0.068mo,  Eq.  (91) 
yields  a  spectral  density 

S6r/r(f)  =  1 .8- 10-8  /fisj  =  S6^/^(f)  =  S6i/i(f),  (92) 

where  I  is  the  source  to  drain  current,  and  N  the  number  of  carriers  in  the  channel.  The  expected 
rms  drain  current  fluctuation  amplitude  is  thus  1.34*  lO"^  Id/NVf.  Both  through  N  and  through  the 
drain  current  this  simple  result  depends  on  the  applied  gate  voltage  Vq. 

7.  Conclusions 


The  preceding  discussion  indicates  that  a  paraelectric  or  hysteresis-free  ferroelectric  gate 
insulation  is  ideal  for  special  HFETs  meeting  bandwidth  requirements  from  0  to  100  GHz.  It 
allows  total  suppression  of  gate  leakage  currents,  while  also  assuring  a  large  increase  in  the 
transconductance  of  the  device  at  frequencies  under  100  MHz,  where  the  permittivity  is  still  high. 

The  gradual  decrease  of  e  in  the  UHF  region  is  actually  very  useful,  since  it  limits  the  free 

fall  of  the  input  impedance  (jcoC)'^  of  the  device  to  zero,  which  would  load  the  source  excessively. 
Finally,  the  lower  temperature  BaMgF4  technology  avoids  oxidation  and  degradation  of  the 

compound  semiconductors,  and  a  further  improvement  in  e  through  mobile  ferroelectric  micro¬ 
domains  is  possible  by  optimizing  the  growing  conditions. 

The  author  is  indebted  to  F.  Schuermeyer  for  suggesting  the  research  for  elimination  of  gate 
currents  and  for  many  helpful  discussions. 
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IX.  TWO-DIMENSIONAL  ALL-OPTICAL  PARALLEL-TO  SERIES 
TDM  AND  DE-MULTIPLEXING  BY 
SPECTRAL-HOLOGRAPHIC  METHODS 

All-optical  one-dimensional  parallel-to  serial  conversion  by  holographic  one-dimensional 
space-to  time  frequency  encoding  was  recently  introduced  [49],  along  with  the  previously  known 
series-to-parallel  transformation.  However,  the  large  discrepancy  between  the  Gbit/s  capacity  of 
coaxial  cables  and  the  Tbit/s  rates  achievable  in  optical  fibers  for  photonic  networks,  requires  the 
multiplexing  of  roughly  1(P  incoming  conventional  signals  in  a  single  optical  fiber.  This  requires 
harnessing  of  the  full  power  of  Fourier-optical  holographic  methods,  by  using  both  dimensions  of 
optical  wavefronts  for  data  processing.  The  processor  system  we  propose  consists  of  two 
independent  optical  channels  Ct  and  Cs  for  carrying  temporal  signals  and  spatial  information,  both 
in  its  parallel-to-series  (Fig.  1")  and  in  its  series-to-parallel  (Fig.  2")  parts,  Ct  is  shown  with 
continuous  lines,  representing  ultra-short  pulsed  beams,  while  Cs  is  shown  with  dotted  rays,  using 
longer  pulses. 

Ct  passes  the  ultrashort  pulse  wave-fronts  through  a  blazed  2-dimensional  reflecting 
diffraction  grating  Gi  which  is  described  below,  followed  by  a  lens  L]  of  focal  length  F,  a  real¬ 
time  hologram  H  in  its  focal  plane,  an  identical  lens  L2,  and  a  second,  similar,  2-D  reflection 
grating  G2,  with  a  distance  F  between  each  of  these  5  elements  on  their  common  optical  axis. 

Cs  is  an  optical  Fourier  transform  channel  which  inputs  the  modified  image  of  a  rectangular 
array  of  sources.  Cs  introduces  this  image  through  the  semi-transparent  mirror  M  together  with  a 
reference  beam.  They  are  introduced  backwards  through  the  lens  L2,  creating  a  Fourier-transform 
hologram  H  which  in  turn  acts  on  the  forward-propagating  Ct  beam,  splitting  it  up  in  a  number  of 
temporally  displaced  beams  equal  with  the  number  of  sources  present  in  the  input  of  Cs-  The  sum 
of  these  beams  is  Fourier-transformed  back  to  the  time-domain,  yielding  the  desired  serial  output 
of  the  given  parallel  data  input.  The  modification  hikes  line  intervals  n  times. 

First  a  general  analysis  of  the  2-D  parallel-to-series  converter  will  be  presented,  and  then  a 
brief  discussion  of  the  corresponding  receiver.  The  multiplexer  works  in  two  steps. 
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Fig.  3:  Two-dimensional  blazed  diffraction  grating,  front  view  and  perspective. 

1)  In  a  first  step,  starting  from  a  single  mode-locked  laser,  a  2-dimensional  manifold  of 
N=nxm  coherent  beams  is  created  by  multiple  beam-splitting,  along  with  a  reference  beam.  By 
synchronized  sampling  and  birefringent  synchronous  modulation,  we  transfer  the  information  from 

each  one  of  the  N  incoming  coaxial  cables  to  one  of  the  N  beams,  in  the  form  of  pulses  of  length  x. 
The  reference  beam  always  carries  the  standard,  unmodulated  pulse  generated  by  the  laser.  Part  of 
the  reference  beam  is  separated  by  a  beam  splitter  and  is  then  re-shaped  [50,51]  for  Ct  to  be  ultra 

short  of  duration  Tq  («x,  up  to  10"3  x,  of  the  order  of  the  reciprocal  carrier  bandwidth,  a  few  tens 


of  femtoseconds).  The  rest  of  the  reference  beam  separated  by  the  beam  splitter  goes  to  the  Cs 
channel  together  with  the  N  modulated  beams.  The  first  step  thus  creates  a  rectangular  array  of 
luminous  sources  with  n  columns  and  m  rows  plus  the  reference  beam,  which  yield  the  (r)  input 

beam  to  Cs,  and  an  ultrashort  sequence  of  pulses  which  yields  the  (tq)  input  beam  to  Ct-  The 
input  to  Ct  is  thus  a  regular  sequence  of  pulses  in  time,  of  the  form 

s(t)  =  p(t-to)exp(jtOct),  (93) 

with  Fourier  transform 


P(to-a)c)exp[-j(o>-(Oc)to],  (94) 

where  P(a))  is  the  Fourier  transform  of  p(t). 


2)  in  the  second  step,  this  input  beam  passes  through  Ct.  Let  the  temporal  transfer 
function  of  Ct  be  H((o).  Then  the  output  spectrum  will  be 


H(a))P((o-(Oc)exp[-j(a)-a)c)tol.  (95) 

To  calculate  H((u),  we  calculate  the  output  resulting  from  the  passage  of  a  monochromatic 
input  exp[j(tot-kr)  through  Ct.  Let  k  be  incident  on  the  2-D  diffraction  grating  located  in  the  xy 

plane  with  directing  cosines  a  and  p  to  the  x  and  y  axes.  The  blazed  2-D  diffraction  grating  is 
obtained  by  applying  a  squared  lattice  chess  board-like  black  deposit  on  a  1-D  blazed  diffraction 
grating  at  45°  to  the  diffraction  grating,  i.e.,  so  that  the  x  and  y  axes  of  the  applied  chess  board  are 
at  45°  from  the  groves  of  the  1-D  blazed  grating  (Fig.  3).  The  lattice  constant  of  the  applied  chess¬ 
board  lattice  is  twice  the  side  of  any  white  or  black  field  on  the  chess  board,  and  equals  the 
constant  of  the  original  blazed  lattice  times  2^^^.  The  diffracted  optical  field  (Fig.  f)  generated  by 
the  monochromatic  input  right  after  the  plane  Pi  of  the  grating  is 

si(x,y;  (d)  =  exp[-j(to-(Oc)(ax+py)/c]  w(x,y),  (96) 

where  w(x,y)  is  the  pupil  function  of  the  grating.  After  passage  through  the  lens  Li  the  optical 
field  in  plane  P2  right  before  the  hologram  is  the  Fourier  transform  of  the  field  in  Eq.  (96),  given 
by 

s2(|,T);  (o)  =  W[|+a(a)-(0c)/2jtc,  Ti+P(a)-o)c)/2jtc],  (97) 

where  W(|,ti)  is  the  spatial  Fourier  transform  of  w(x,y).  The  wave  numbers  |  and  ti  are 
expressed  in  terms  of  the  spatial  Cartesian  coordinates  X  and  Y  in  P2  by  §=oi)X/2jicF  and 

Ti=a)Y/2jtcF.  We  note  parenthetically  that,  as  indicated  by  Eq.  (97),  an  optical  pulse  introduced 
into  the  system  will  have  its  various  Fourier  components  spatially  dispersed  along  the  diagonal  of 
the  Fourier  transform  plane,  each  occupying  an  area  determined  by  the  function  W(§,ti)  . 

We  now  take  into  account  that  a  spatial  Fourier  transform  hologram  of  the  above-mentioned 
rectangular  nxm  input  array  has  been  independently  recorded  by  Cs  in  the  hologram  H  of  Fig.  l'! 
This  hologram  is  a  superposition  of  the  nxm  holograms  of  the  individual  luminous  sources  present 
in  the  rectangular  input  array,  with  amplitudes  Afxv  corresponding  to  the  bit  of  information  carried 

by  the  channel  |xv  to  be  multiplexed.  This  hologram  serves  as  a  temporal  frequency  filter  with 
transmittance 
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(98) 


m  n 

t(|,Ti)=  2  2  AjA,vexp|j2jt(cow/w)(vAxi+^AyT])], 

p=lv=l 

where  Ax  is  the  geometrical  distance  between  columns  and  Ay  is  the  modified  distance  between  the 
rows  in  the  regular  input  array  fed  into  Cs,  whereas  (Ow  is  the  frequency  of  the  writing  field  used 
for  recording  the  hologram,  in  our  case  coc-  It  is  necessary  to  take  Ay  ^  nAx,  and  suitable  to 

choose,  e.g,.  Ay  =  (n+l)Ax.  This  is  achieved  with  the  help  of  a  convergent  cylindrical 
magnifying  lens  placed  in  front  of  the  rectangular  array  of  light  sources,  with  its  cylindrical  axis 
parallel  to  the  x-axis,  which  would  cause  the  separation  between  rows  to  appear  larger.  In 

addition,  we  assume  (n+l)mAx+cTo  ^  ex. 

Returning  now  to  Q,  right  after  passage  through  the  hologram,  the  optical  field  in  P3  will 
be 

ssd,!!;  co)  =  s2(§,'n;  to)  t(|,ri) .  (99) 

After  passage  through  L2,  the  field  in  plane  P4  is  given  by  a  second  spatial  Fourier  transform 
which  yields 

m  n 

S4(x,y;  to)  =  ^  2  A^,vw[-x+vAx,  -y+pAyJ 
^=lV=:l 

xexp[(j/c)(to-(Oc)[a(x-vAx(o)w/(o))  +  P(y-pAy((Ow/to))].  (100) 

This  field  is  diffracted  by  the  2-D  grating  G2  and  yields  the  output  field  propagating  in  the 
k"  or  z"  direction 

m  n 

s5(x",y”;  to)  =  ^  2  Ajx,v'''[-x"+ct"vAx,  -y"+P"pAy] 

p^=lv=l 

xexp[(j/c)((o-(Oc)(tOw/(o)(avAx  +  ppAy)]  =  H(o>).  (101) 

This  is  the  temporal  transfer  function  of  the  system.  Here  a"  and  (3"  are  directing  cosines  of  the 

outgoing  wave  vector  k".  Assuming  Ay  =  (n+l)A  and  A=Ax, 
m  n 

H(to)  =22  A[x,vw[-x"+a"vAx,  -y"+p>Ay] 
p=lv=l 

xexp{j(A/c)((o-tOc)(tOw/(o)[av  +  Pp(n+1)]}  (102) 

The  temporal  output  function  allows  the  determination  of  the  serial  output  of  the  multiplexer 
as  the  Fourier  transform  of  Eq.  (95) 
m  n 

So(x",y";  t)  =  [  2  2  AjA,vw[-x"+a"vAx,  -y"+p>Ay]p(t-fo  -r6t)]  expCjtOcO.  (103) 

(X=lv=l 

Here  r  =  v  +  fx(pAy/aAx)  =  v  +  p.(p/a)(n+l);  for  a=p,  r  =  v  +  p(n+l).  The  interval  6t  was 

defined  as  <xA\/c.  The  ratio  cowAo  was  approximated  by  1  which  is  reasonable  for  not  too  large 
bandwidths,  but  may  lead  to  errors  in  our  case.  These  errors  should  be  further  investigated. 

Eq.  (103)  shows  that  the  parallel-to-series  or  multiplexing  operation  was  performed.  Along 
with  the  series  signal,  the  non-diffracted  portion  of  the  readout  beam  coming  from  the  hologram  H 
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in  Fig.  1  is  also  transmitted  to  the  receiver  as  a  reference  signal. 

In  principle,  all  beams  could  be  reverse^.  This  observation  makes  the  description  and 
operation  of  the  de-multiplexer  shown  in  Fig.  2:  OTaightforward.  A  hologram  is  recorded  at  the 
receiver  site  between  the  incoming  serial  signal  shown  by  Eq.  (103)  and  the  received  reference 
signal,  both  having  first  been  diffracted  in  reflection  by  the  identical  diffraction  grating  G.  This 
real  time  hologram  is  modulating  the  cw  monochromatic  readout  wave  front  coming  from  the  left  in 

if  ^  ® 

Fig.  2  and  passing  through  the  hologram,  thereby  creating  N=nxm  beams.  The  result  is  thus  a 
two-dimensional  nxm  parallel  optical  output  which  is  extracted  with  the  help  of  the  semi¬ 
transparent  mirror  M  and  projected  on  a  nxm  array  of  outgoing  optical  fibers,  or  on  a  similar  array 
of  photodetectors. 
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